AFIT/GA/ENY/93D-5 


DTIC 

electe 

Pi 


i 


AD-A273  723 

lllllllll 


AN  APPROXIMATE  SOLUTION  FOR  THE 
SPINUP  DYNAMICS  OF  NEAR  AXISYMMETRIC  AXIAL 
GYROSTATS  USING  THE  METHOD  OF  MULTIPLE  SCALES 

THESIS 

Stewart  J  Kowall 
Captain,  USAF 

AFIT/GA/ENY/93D-5 


93-30473 


Approved  for  public  release;  distribution  unlimited 

93  12  15  083 


AFIT/GA/ENY/93D-5 


AN  APPROXIMATE  SOLUTION  FOR  THE  SPINUP  DYNAMICS 
OF  NEAR  AXISYMMETRIC  AXIAL  GYROSTATS  USING  THE 
METHOD  OF  MULTIPLE  SCALES 


THESIS 


Presented  to  the  Faculty  of  the  Graduate  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 


Air  University 
In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science  in  Astronautical  Engineering 


Stewart  J  Kowall,  B.S.  Aeronautical  Engineering 
Captain,  USAF 


December,  1993 


Accesion  For 

NTIS  CRA&I 
DTIC  TAB 
Unannounced 
Justification 


□ 

□ 


By . . . 

Distribution  f 


Availability  Codes 


Oist 


Avail  and/or 
Special 


0  ^ 

lie**! 


It 


Approved  for  public  release;  distribution  unlimited 


Acknowledgements 

The  work  contained  in  this  thesis  is  the  result  of  mzuiy  hours  of  work  which 
would  not  have  been  completed  but  for  the  help  of  three  groups  of  people.  The  first  is 
my  advisor,  Capt  Christopher  D.  Hall.  Without  his  patience  and  understanding,  this 
thesis  would  ,  most  likely,  not  have  been  completed  on  time.  Capt  Hall  is  an  excellent 
instructor.  Should  he  ever  decide  to  leave  his  teaching  position,  it  would  be  a  great 
loss  for  his  future  students,  but  a  wonderful  gain  for  his  next  employer.  Second,  I 
would  like  to  thank  Lt  Col  William  P.  Baker.  His  long  hours  of  extra  instruction  and 
contribution  to  this  thesis  taught  me  a  great  deal  about  the  nature  of  mathematics. 
Finally,  I  say  “Thank  You”  to  my  family,  Janet,  Nathan  and  Nicholas.  The  price 
paid  by  you  in  terms  of  hours  spent  in  the  books  versus  hours  spent  at  home  was 
very  large.  I  only  hope  that  I  can  repay  you  in  the  years  to  come. 

Stewart  J  Kowall 


II 


Table  of  Contents 


Page 


Acknowledgements .  ii 

List  of  Figures  .  v 

List  of  Tables .  x 

Abstrect .  xi 

I.  Introduction .  1-1 

1.1  Problem  Statement .  1-3 

1.2  Related  Work .  1-3 

1.3  Outline  of  Thesis .  1-5 

II.  Dynamics  of  the  Dual-Spin  Spacecraft  .  2-1 

2.1  Dimensionless  Equations .  2-5 

2.2  SpinUp  Bo\mdary  Conditions .  2-6 

2.3  The  Dynamical  Shape  of  the  Gyrostat .  2-7 

2.4  Spinup  Problems .  2-10 

2.5  Summary .  2-14 

III.  The  Straightforward  Expansion .  3-1 

3.1  Straightforward  Expansion  Results .  3-8 

3.2  Conclusions .  3-13 

IV.  The  Method  of  Multiple  Scales  .  4-1 

4.1  Method  of  Multiple  Scales  Results .  4-15 

4.1.1  Oblate  Gyrostats .  4-16 


111 


Page 

4.1.2  Prolate  Gyrostats .  4-27 

4.1.3  Transverse  Spinup  .  4-31 

4.2  Summary .  4-34 

V.  Cylindrical  Coordinates .  5-1 

5.1  A  Comparison  Between  the  C2irte8ian  and  Cylindrical 

Approximate  Solutions .  5-11 

5.1.1  Oblate  Gyrostats .  5-11 

5.1.2  Prolate  Gyrostats .  5-15 

5.1.3  Transverse  Spinup  .  5-16 

5.2  Summary .  5-20 

VI.  Conclusions  and  Recommendations .  6-1 

6.1  Conclusions .  6-1 

6.2  Recommendations .  6-2 

Appendix  A.  Straightforward  Expansion  Method  Solution  Terms  .  .  A-1 

Appendix  B.  FVesnel  Function  Integrals .  B-1 

Bibliography  .  BIB-1 

Vita .  VITA-1 


IV 


List  of  Figures 

Figure  Page 

2.1.  Model  of  axial  gyrostat.  The  asymmetric  platform  is  denoted  as  V 
while  the  axisymmetric  rotor  is  denoted  as  H.  The  body  frame 

is  aligned  with  the  principle  axis  and  rotates  relative  to  the  inertial 
frame  Fi.  The  total  angular  momentum  vector  is  denoted  as  h  and 
is  constant  in  direction  and  magnitude .  2-2 

2.2.  Allowable  values  of  and  13  for  a  gyrostat  with  a  rod  shaped  rotor 

(ii  =  0).  Oblate  gyrostats  have  negative  values  of  12  and  13.  Prolate 
gyrostats  have  positive  values  of  tj  and  13 .  2-8 

2.3.  Momentum  spheres  for  a  near  axisymmetric  gyrostat  with  a  one 

percent  difference  in  transverse  inertias  zmd  a  range  of  /<’s .  2-12 

3.1.  Fresnel  Sine  and  Fresnel  Cosine  functions .  3-6 

3.2.  Straightforward  expansion  method  approximate  solution  of  the  xi 

component  of  angular  momentum  of  an  oblate  gyrostat  with  initial 
conditions  Xi  =  0.9,2]  =  23  =  0.30822  and  e  =  0.01;  gf  =  0.1  .  .  .  3-10 

3.3.  Straightforward  expansion  method  approximate  solution  of  the  xi 

component  of  angular  momentum  of  an  oblate  gyrostat  with  initial 
conditions  21  =  0.9,2]  =  23  =  0.30822  and  e  =  0.01;  g  =  0.01  .  .  3-10 

3.4.  Straightforward  expansion  method  approximate  solution  of  the  X] 

component  of  angular  momentum  of  an  oblate  gyrostat  with  initial 
conditions  21  =  0.9,2]  =  23  =  0.30822  and  e  =  0.01;  =  0.1  .  .  .  3-11 

3.5.  Strsdghtforward  expansion  method  approximate  solution  of  the  X] 

component  of  angular  momentum  of  an  oblate  gyrostat  with  initial 
conditions  2]  =  0.9,2]  =  23  =  0.30822  and  e  =  0.01;  g  =  0.01  .  .  3-11 

3.6.  Maximum  error  detected  in  the  Xj  component  of  angular  momen¬ 

tum  as  a  function  of  t/inai  (or  g)  for  an  oblate  gyrostat  using  the 
straightforward  expansion  method .  3-12 

V 


Page 


Figure 

3.7.  Maximum  error  detected  in  the  Xj,  X3  component  of  angular  mo¬ 
mentum  as  a  function  of  t final  (or  g)  for  an  oblate  gyrostat  using 
the  straightforward  expansion  method .  3-12 

4.1.  zii  Variation  for  an  oblate  gyrostat  with  initi2d  conditions  Xi  = 

0.7,12  =  0.4, 13  =  0.59161  and  e  =  0.01;  g  =  0.01 .  4-11 

4.2.  Numerical  vs  approximate  solution  of  the  Xi  component  of  angu¬ 
lar  momentum  of  an  oblate  gyrostat  with  initial  conditions  xi  = 

0.7,X2  =  0.4,  X3  —  0.59161  and  e  =  g  =  0.01 .  4-17 

4.3.  Percent  relative  error  of  the  Xi  component  of  angular  momentum 
of  an  oblate  gyrostat  with  initial  conditions  xi  =  0.7, X2  =  0.4, 

X3  =  0.59161  and  e  =  g  =  0.01  .  4-17 

4.4.  Numerical  vs  approximate  solution  of  the  X2  component  of  angu¬ 
lar  momentum  of  an  oblate  gyrostat  with  initial  conditions  xi  = 

0.7,X2  =  0.4,  X3  =  0.59161  and  £  —  g  =  0.01 .  4-18 

4.5.  Percent  relative  error  of  the  X2  component  of  angular  momentum 
of  an  oblate  gyrostat  witn  initial  conditions  xi  =  0.7,X2  =  0.4, 

X3  =  0.59161  and  e  =  g  =  0.01  .  4-18 

4.6.  Numerical  vs  approximate  solution  of  the  X3  component  of  angu¬ 
lar  momentum  of  an  oblate  gyrostat  with  initial  conditions  Xi  = 

0.7,X2  =  0.4,  X3  =  0.59161  and  e  =z  g  =  0.01 .  4-19 

4.7.  Percent  relative  error  of  the  X3  component  of  angular  momentum 
of  an  oblate  gyrostat  with  initial  conditions  xi  =  0.7,X2  =  0.4, 

X3  =  0.59161  and  £  —  g  =  0.01  .  4-19 

4.8.  Amplitude  error  for  the  X2  component  of  angular  momentum  of  an 
oblate  gyrostat  with  initial  conditions  xi  =  0.7,  X2  =  0.4,  X3  = 

0.59161  and  £  =  if  =  0.01 .  4-21 

4.9.  xii  For  two  diflferent  initial  conditions .  4-21 

4.10.  Maximum  error  detected  as  a  function  of  initial  conditions  for  sm 

oblate  gyrostat .  4-22 

4.11.  Maximum  error  detected  in  the  last  10  percent  of  spinup  as  a  func¬ 
tion  of  initial  conditions  for  an  oblate  gyrostat .  4-23 


VI 


Page 


Figure 

4.12.  Numerical  vs  approximate  solution  of  the  X2  component  of  angu¬ 
lar  momentum  of  an  oblate  gyrostat  with  initial  conditions  Xi  = 

0.7, xj  =  0.4,  X3  =  0.59161  and  e  =  0.01;^  =  0.001 .  4-24 

4.13.  Numerical  vs  approximate  solution  of  the  X3  component  of  2mgu- 
lar  momentum  of  an  oblate  gyrostat  with  initial  conditions  Xi  = 

0.7,X2  =  0.4,  X3  =  0.59161  and  e  =  0.01;  =  0.001 .  4-25 

4.14.  Phase  error  i,hat  develops  for  an  oblate  gyrostat  with  g  «e.  .  .  4-25 

4.15.  Validity  of  solutions  with  g  «e .  4-26 

4.16.  Maximum  error  detected  as  a  function  of  c  for  an  oblate  gyrostat 

with  initial  conditions  x\  =  0.7,X2  =  0.4,  X3  =  0.59161  ;g  =  e  .  .  .  4-27 

4.17.  Cartoon  of  difiFerent  trajectory  paths  in  a  ^  =  0  nonlinear  system.  4-29 

4.18.  Numerical  vs  approximate  solution  of  the  Xi  component  of  angu¬ 
lar  momentum  of  a  prolate  gyrostat  with  initial  conditions  Xi  = 

0.9,X2  =  0.30822,  X3  =  0.30822  and  e  =  5  =  0.01 .  4-30 

4.19.  Numerical  vs  approximate  solution  of  the  Xj  component  of  angu¬ 
lar  momentum  of  a  prolate  gyrostat  with  initial  conditions  Xi  = 

0.9, X2  =  0.30822,  X3  =  0.30822  and  e  =  g  =  0.01 .  4-30 

4.20.  Numerical  vs  approximate  solution  of  the  X3  component  of  angu¬ 
lar  momentum  of  a  prolate  gyrostat  with  initial  conditions  Xi  = 

0.9,X2  =  0.30822,  X3  =  0.30822  and  e  =  g  =  0.01 .  4-32 

4.21.  Percent  relative  error  in  the  Xi  component  of  angular  momentum 
of  a  prolate  gyrostat.  Initial  conditions  are  matched  at  ^  =  0  and 

H  =  0.505 .  4-32 

4.22.  Percent  relative  error  in  the  X2  component  of  angular  momentum 
of  a  prolate  gyrostat.  Initial  conditions  are  matched  at  ^  =  0  and 

fi  =  0.505 .  4-33 

4.23.  Percent  relative  error  in  the  X3  component  of  angular  momentum 
of  a  prolate  gyrostat.  Initial  conditions  are  matched  at  /x  =  0  and 

fi  =  0.505 .  4-33 

4.24.  Numerical  vs  approximate  solution  of  the  xi  component  of  angular 

momentum  of  a  prolate  gyrostat  during  transverse  spinup .  4-34 


Vll 


I- 


Figure  Page 

4.25.  Numerical  vs  approximate  solution  of  the  X]  component  of  angular 

momentum  of  a  prolate  gyrostat  during  transverse  spinup .  4-35 

4.26.  Numerical  vs  approximate  solution  of  the  X3  component  of  angular 

momentum  of  a  prolate  gyrostat  during  transverse  spinup .  4-35 

4.27.  Numerical  vs  approximate  solution  of  the  xx  component  of  angular 

momentum  of  a  prolate  gyrostat  during  transverse  spinup .  4-36 

4.28.  Numerical  vs  approximate  solution  of  the  X}  component  of  angular 

momentum  of  a  prolate  gyrostat  during  transverse  spinup.  ....  4-36 

4.29.  Numerical  vs  approximate  solution  of  the  X3  component  of  angular 

momentum  of  a  prolate  gyrostat  during  transverse  spinup.  ....  4-37 


5.1.  Comparison  between  cartesian  and  cylindrical  equations  for  per¬ 
cent  relative  error  of  the  Xi  component  of  angular  momentum  of 


an  oblate  gyrostat  with  initial  conditions  Xx  =  0.7,  X3  =  0.4,  X3  = 

0.59161  and  e  =  g  =  0.01 .  5-13 


5.2.  Comparison  between  cartesian  and  cylindrical  equations  for  per¬ 
cent  relative  error  of  the  xj  component  of  angular  momentum  of 
an  oblate  gyrostat  with  initial  conditions  Xx  =  0.7,  Xj  =  0.4,  X3  = 

0.59161  and  e  =  g  =  0.01 .  5-13 

5.3.  Comparison  between  Ccirtesian  and  cylindrical  equations  for  per¬ 
cent  relative  error  of  the  X3  component  of  angular  momentum  of 
an  oblate  gyrostat  with  initial  conditions  xx  =  0.7,  X2  =  0.4,  X3  = 

0.59161  and  e  =  g  =  0.01 .  5-14 

5.4.  Comparison  between  cartesian  and  cylindrical  equations  for  per¬ 
cent  relative  error  of  the  xj  component  of  angidar  momentum  of  an 
oblate  gyrostat  with  initial  conditions  xx  =  0.4,  xj  =  0.64807,  X3  = 

0.64807  and  e  =  g  =  0.01 .  5-14 

5.5.  Comparison  between  cartesian  and  cylindrical  equations  for  per¬ 
cent  relative  error  of  the  X3  component  of  angular  momentum  of  am 
oblate  gyrostat  with  initial  conditions  xx  =  0.4,  X2  =  0.64807,  X3  = 

0.64807  and  e  =  g  =  0.01 .  5-15 


vm 


Figure  Page 

5.6.  Comparison  between  cartesian  and  cylindrical  equations  of  maxi¬ 
mum  error  detected  in  the  last  ten  percent  of  spin  up  as  a  function 

of  initial  conditions  for  an  oblate  gyrostat .  5-16 

5.7.  Numerical  vs  cylindrical  and  cartesian  coordinate  approximate  so¬ 
lution  of  the  Xi  component  of  angular  momentum  of  a  prolate  gy¬ 
rostat  with  initial  conditions  xi  =  0.4,  xj  =  0.64807,  X3  =  0.64807 

and  e  =  =  0.01 .  5-17 

5.8.  Numerical  vs  cylindrical  and  cartesian  coordinate  approximate  so¬ 
lution  of  the  X}  component  of  angular  momentum  of  a  prolate  gy¬ 
rostat  with  initial  conditions  Xi  =  0.4,  xj  =  0.64807,  X3  =  0.64807 

and  e  =  g  =  0.01 .  5-17 

5.9.  Numerical  vs  cylindrical  and  cartesian  coordinate  approximate  so¬ 
lution  of  the  X3  component  of  angular  momentum  of  a  prolate  gy¬ 
rostat  with  initial  conditions  Xi  =  0.4,  X2  =  0.64807,  X3  =  0.64807 

and  e  =  g  =  0.01 .  5-18 

5.10.  Percent  relative  error  of  the  Xi  component  of  angular  momentum 

of  a  prolate  gyrostat  using  cylindrical  coordinates  with  initial  con¬ 
ditions  xi  =  0.7,  X2  =  0.4,  *3  =  0.59161  and  e  =  g  ~  0.01 .  5-18 

5.11.  Percent  relative  error  of  the  xj  component  of  angular  momentum 

of  a  prolate  gyrostat  using  cylindrical  coordinates  with  initial  con¬ 
ditions  x\  =  0.7, 12  =  0.4,  X3  =  0.59161  and  e  =  g  =  0.01 .  5-19 

5.12.  Percent  relative  error  of  the  X3  component  of  angular  momentum 

of  a  prolate  gyrostat  using  cylindrical  coordinates  with  initial  con¬ 
ditions  xi  =  0.7,  xj  =  0.4,  X3  =  0.59161  and  e  =  g  =  0.01 .  5-19 

5.13.  Numerical  vs  cylindrical  and  cartesian  coordinate  approximate  solu¬ 

tion  of  the  xi  component  of  angular  momentum  of  a  prolate  gyrostat 
during  transverse  spinup .  5-20 

5.14.  Numerical  vs  cylindrical  and  cartesian  coordinate  approximate  solu¬ 

tion  of  the  xj  component  of  angular  momentum  of  a  prolate  gyrostat 
during  transverse  spinup .  5-21 

5.15.  Numerical  vs  cylindrical  and  csurtesian  coordinate  approximate  solu¬ 

tion  of  the  X3  component  of  angular  momentum  of  a  prolate  gyrostat 
during  transverse  spinup .  5-21 


IX 


List  of  Tables 


Table  Page 

2.1.  Equilibrium  Points .  2-11 


r 


AFIT/GA/ENY/93D-5 

Abstract 

Approximate  solutions  for  the  spinup  of  a  near  axisymmetric  dual-spin  gy¬ 
rostat  are  derived  using  the  straightforward  expansion  method  and  the  method  of 
multiple  scades.  Two  method  of  multiple  scale  solutions  are  presented.  The  first 
is  derived  using  cartesian  coordinates  while  the  second  is  derived  using  cylindrical 
coordinates.  The  multiple  scales  solutions  are  compared  to  numerically  integrated 
results  for  oblate  and  prolate  configurations.  A  comparison  for  flat  spin  recovery  is 
also  accomplished.  Excellent  results  are  obtained  for  oblate  configurations.  Trajec¬ 
tory  separatrix  crossings  hindered  the  results  for  prolate  configurations  and  flat  spin 
recoveries. 
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AN  APPROXIMATE  SOLUTION  FOR  THE  SPINUP  DYNAMICS 
OF  NEAR  AXISYMMETRIC  AXIAL  GYROSTATS  USING  THE 
METHOD  OF  MULTIPLE  SCALES 


1.  Introduction 

In  1958,  the  United  States  of  America  launched  Explorer  I  into  orbit  as  its 
answer  to  the  Soviet  Unions’  Sputnik.  The  great  space  race  between  the  two  nations 
began,  and  unfortunately,  the  United  States’  entry  suffered  a  serious  flaw.  Explorer 
I  was  shaped  like  a  missile  and  was  designed  to  rotate  about  its  minimum  moment 
of  inertia  axis.  Ninety  minutes  after  insertion  into  its  orbit.  Explorer  I  was  tumbling 
end-over-end  and  the  engineers  at  the  Jet  Propulsion  Laboratory  did  not  know  why 
(9:129).  In  the  following  years,  U.S.  engineers  learned  a  great  deal  about  rotational 
dynamics. 

By  1960,  engineers  accepted  the  idea  that  internal  energy  dissipation,  caused 
by  spacecraft  flexibility,  violated  the  assumptions  of  the  classical  analyses  of  Euler. 
The  result  was  the  belief  that  all  spin  stabilized  satellites  must  be  shaped  like  “tuna 
cans”  so  that  spin  stabilization  could  be  achieved  through  spin  about  the  maximum 
moment  of  inertia  axis.  As  satellite  design  progressed,  the  idea  of  an  inertial  platform 
became  increasingly  attractive.  With  such  a  platform,  sensors  could  be  continuously 
pointed  at  a  source  rather  than  rely  on  scanned  coverage.  The  first  satellite  to  employ 
such  a  platform  was  launched  in  1962  and  was  called  the  Orbiting  Solar  Observatory, 
080-1(13:761).  This  satellite  was  still  designed  to  spin  about  its  maximum  moment 
of  inertia  axis,  but  the  concept  of  an  inertial  platform  was  flown  (9:131). 

As  satellites  progressed  even  further,  however,  they  were  designed  to  accom¬ 
plish  more  tasks  and  as  a  result,  became  bigger.  The  constraint  to  such  large  satellites 
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was  the  size  of  the  payload  shroud  used  to  protect  them  during  assent  to  orbit.  The 
aerodynamics  of  rocket  design  dictated  a  long  slender  shape  which  conflicted  with 
the  short  squat  shape  required  for  spin  stability.  The  lessons  of  Explorer  I,  however, 
were  firmly  entrenched  in  the  minds  of  our  satellite  designers.  Surprisingly,  Mr. 
Vernon  Landon,  an  employee  at  RCA,  deduced  the  instability  about  the  minimum 
moment  of  inertia  axis  due  to  energy  dissipation  as  early  as  1957.  Further,  by  1962, 
Mr.  L2uidon  knew  that  stable  spin  about  the  minimum  moment  of  inertia  axis  was 
possible  if  enough  angular  momentum  was  stored  in  a  rigid  rotor  aligned  with  the 
spin  axis  of  the  spacecraft  (9:130).  Unfortunately,  his  efforts  went  unrecognized  for 
mamy  years. 

In  1965,  Mr.  Tony  lorillo  of  Hughes  Aircraft  Company,  found  that  by  putting 
a  damper  on  the  platform  section  of  a  dual-spin  spacecraft,  spin  stabilization  for 
spacecraft  of  any  inertia  distribution  was  possible.  By  using  an  energy  sink  amalysis, 
it  was  shown  that  if  the  energy  dissipation  rate  of  the  platform  is  much  greater  than 
the  energy  dissipation  rate  of  the  rotor,  spin  stabilization  could  be  achieved  for  a 
spacecraft  spinning  about  its  minor  moment  of  inertia  axis  (14:151-15^). 

The  advantages  of  dual-spin  spacecraft  are  many  and  it  is  used  on  many  of 
todays  satellites.  There  is  an  important  problem  in  the  dynamics  of  dual-spin  space¬ 
craft  to  which  this  thesis  addresses  itself.  That  problem  is  concerned  with  the  spinup 
of  the  rotor  such  that  spin  stability  can  be  achieved.  The  spinup  memeuver  can  occur 
at  two  different  instances  in  the  spacecrafts  lifespan.  The  first  is  when  a  satellite  is 
initially  deployed.  As  the  satellite  is  inserted  into  its  orbit,  the  rotor  and  platform 
are  both  spinning  at  the  same  rate  and  their  relative  angular  velocity  is  zero.  In  order 
to  give  the  platf  rm  an  essentially  inertial  reference,  a  torque  motor  is  incorporated 
to  spinup  the  rotor  amd  thereafter  maintain  the  desired  relative  rotation  rate.  The 
initicd  spinup  causes  the  angular  momentum  of  the  platform  to  be  transferred  to 
the  rotor  which,  in  turn,  causes  the  platform  to  despin.  The  second  instance  occurs 
if  the  torque  motor  used  to  maintain  a  specified  spin  rate  between  the  rotor  and 
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platform  fails.  Internal  friction  will  eventually  reduce  the  relative  spin  rate  between 
the  rotor  and  platform  to  zero.  Continued  energy  dissipation  will  cause  the  satellite 
to  seek  its  maximum  moment  of  inertia  axis.  Should  failure  of  the  torque  motor  be 
corrected,  the  relative  spin  rate  between  the  rotor  and  platform  and  the  satellite’s 
attitude  could  be  reestablished  by  spinning  up  the  rotor. 

1.1  Problem  Statement 

The  objective  of  this  thesis  is  to  derive  an  approximate  solution  to  the  equa¬ 
tions  of  motion  that  describe  the  spinup  dynamics  of  near  axisymmetric  dual-spin 
satellites.  Specifically,  the  method  of  multiple  scciles  perturbation  technique  will  be 
used  to  develop  these  approximate  solutions  valid  for  dual-spin  satellites  with  small 
differences  in  transverse  inertias  being  spun  up  with  a  non-zero  axial  torque.  There 
are  several  beneffts  associated  with  the  development  of  approximate  solutions.  The 
first  is  the  reduction  of  computational  time.  If  the  attitude  of  a  spacecraft  is  desired 
for  a  particular  time  in  the  future,  the  equations  of  motion  must  be  integrated  from 
the  initial  state  to  the  time  in  question.  If  an  approximate  solution  is  available,  the 
time  need  only  be  placed  in  the  equation  to  determine  the  future  attitude.  Second, 
approximate  solutions  provide  valuable  insight  into  the  dynamics  of  a  rotating  body 
as  its  trajectory  passes  through  a  separatrix.  Third,  with  an  increased  level  of  under¬ 
standing  of  the  spinup  dynamics  attained,  weight  s  wings  for  future  spacecraft  may 
be  resJized  since  torque  values  for  spinup  motors  can  be  more  accurately  predicted. 
Finally,  development  of  approximate  solution  techniques  further  enhance  the  body 
of  knowledge  available  for  the  design  of  future  space  systems  (3:860). 

1.2  Related  Work 

Four  different  groups  of  researchers  have  completed  work  in  approximate  so¬ 
lutions  that  is  most  closely  related  to  the  type  of  dual-spin  satellite  modelled  in 
this  thesis.  Sen  and  Bainum  (13)  derived  an  approximate  solution  for  the  spinup 
of  an  oblate  satellite.  Using  a  perturbation  technique,  they  developed  dimensional 
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equations  for  a  near  axisymmetric  gyrostat  model  of  the  Small  Astronomy  Satellite 
(SAS-A).  Their  results  only  included  the  transverse  components  of  angular  momen¬ 
tum.  While  they  were  able  to  capture  the  amplitude  of  the  system,  they  were  not 
able  to  account  for  the  increase  in  frequency  during  the  spinup  maneuver. 

The  work  of  Gebman  and  Mingori  (3)  was  completed  in  1976  and  was  concerned 
with  an  approximate  solution  for  a  flat  spin  recovery.  Their  solution  incorporated 
the  use  of  a  multiple  scales  method  in  which  their  initial  conditions  were  fixed  at  the 
exact  equilibrium  point  for  flat  spin  recovery.  Their  solution  proved  satisfactory  for 
this  one  particular  case. 

A  method  by  which  the  spinup  dynamics  of  a  dual-spin  satellite  in  an  external 
torque  free  environment  could  be  analyzed  has  recently  been  developed  by  Hall  (5). 
By  plotting  the  energy  of  the  satellite  versus  the  angular  momentum  of  the  rotor,  he 
was  able  to  capture  the  dynamics  of  three  distinct  spinup  problems.  The  first  two 
problems  were  oblate  and  prolate  spinup  which  are  associated  with  the  initial  satellite 
deployment.  The  terms  “oblate”  and  “prolate”  refer  to  the  inertia  distribution  of 
the  satellite.  The  third  problem  was  termed  transverse  spinup  and  is  associated  with 
the  satellite  attitude  recovery  problem  should  the  torque  motor  fail.  This  problem  is 
commonly  referred  to  as  flat  spin  recovery.  Hall  took  the  four  first  order  differential 
equations  that  describe  the  dynamics  of  a  rigid  dual-spin  satellite  and  reduced  them 
to  a  single  first-order  non-autonomous  ordinary  differential  equation  for  the  slow 
evolution  of  kinetic  energy  during  spinup.  This  reduction  was  made  possible  by 
using  the  conservation  of  angular  momentum  and  the  method  of  averaging.  Previous 
studies  of  spin-up  dynamics  had  been  restricted  to  particul2u:  classes  of  gyrostats  or 
particular  regions  of  phase  space.  Hall,  however,  provides  a  unified  treatment  for 
all  gyrostats  with  his  discovery  of  a  symmetry  relating  oblate  and  prolate  gyrostats 
which  reduce  the  number  of  cases  to  be  investigated.  A  significant  difference  in 
notation  between  this  thesis  and  Hall’s  work  is  in  the  use  of  the  variable  “e”.  Hall 
used  “e”  to  denote  the  non-dimensional  torque  applied  by  the  platform  on  the  rotor. 
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This  thesis  used  the  variable  “jr”  to  denote  this  value  while  “‘e”  is  used  as  the  small 
parameter  in  the  perturbation  methods. 

Kinney  (8)  linearized  the  equations  of  motion  for  the  spinup  of  an  ideal  ax¬ 
ial  gyrostat  about  the  oblate  and  prolate  equilibrium  points.  He  then  applied  the 
WKB  (Wentzel,  Krammers,  and  Brillouin)  approximation  method  for  the  linearized 
equations  of  motion.  Because  the  equations  of  motion  were  linearized  about  the 
equilibrium  points,  his  approximate  solutions  were  forced  to  stay  very  close  to  these 
points.  While  obtaining  good  results  for  oblate  spinup,  the  WKB  solution  failed  for 
prolate  spinup  due  to  its  inherent  inability  to  2Mlapt  to  the  bifurcations  associated 
with  the  prolate  equilibrium  point. 

l.S  Outline  of  Thesis 

The  thesis  begins  with  a  review  of  the  dynamics  of  a  dual-spin  spacecraft. 
Chapter  2  first  develops  the  equations  of  motion  for  the  spinup  maneuver  and  then 
transforms  them  to  a  non-dimensional  system  of  equations.  The  boundary  conditions 
for  the  spinup  maneuver  are  discussed  and  a  numerical  value  for  the  end  of  spinup 
is  derived.  In  order  to  provide  a  treatment  for  all  types  of  dual-spin  spacecraft,  the 
dynamical  shape  of  the  spacecraft  in  terms  of  its  non-dimensional  moments  of  inertia 
is  discussed.  Finally,  a  review  of  the  three  types  of  spinup  problems  (oblate,  prolate, 
and  transverse)  is  presented. 

Chapter  3  employs  the  straightforward  expansion  method  to  the  equations  of 
motion.  While  this  method  fails  in  general,  it  provides  insight  for  the  development  of 
an  approximate  solution  using  the  method  of  multiple  scales.  Results  for  an  oblate 
spacecraft  are  presented  to  show  the  inadequacy  of  this  expansion. 

Chapter  4  employs  the  method  of  multiple  scales  to  derive  another  approximate 
solution.  Excellent  results  are  obtained  for  oblate  spacecraft.  Results  for  prolate 
spacecraft  and  transverse  spinup  are  also  presented  and  discussed. 


Chapter  5  also  employs  the  method  of  multiple  scales,  but,  applies  the  method 
to  the  equations  of  motion  after  they  are  transformed  to  cylindrical  coordinates. 
Comparison  between  the  two  multiple  scales  approximate  solutions  2ure  presented  for 
the  three  spinup  problems. 
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11.  Dynamics  of  the  Dual-Spin  Spacecraft 

The  spacecraft  to  be  modelled  is  as  shown  in  Fig  2.1  (5:26).  The  model  con¬ 
sists  of  a  rigid  body  free  to  rotate  in  space.  The  relative  motion  of  the  platform 
('P)  and  rotor  (P),  however,  is  constrained  to  rotation  about  a  rigid  frictionless 
shaft  aligned  with  a  principal  axis  (ej).  This  model  accommodates  an  asymmetric 
spacecraft  through  the  assumption  that  the  platform  is  asymmetric,  but  the  rotor 
is  axisymmetric  about  a  principal  axis  which  is  also  the  relative  rotation  axis.  This 
assumption  leads  to  a  constsmt  moment  of  inertia  tensor  for  the  modelled  spacecraft, 
and  therefore,  by  definition  is  also  known  as  a  gyrostat.  The  vector  h  represents  the 
total  angular  momentum  of  the  gyrostat.  Because  the  assumption  of  no  external 
torques  is  made,  the  magnitude  of  the  angular  momentum  vector  must  be  conserved 
and  the  angular  momentum  vector  is  fixed  in  an  inertial  reference  frame. 

=  hl  +  hl  +  hl  =  constant  (2.1) 

In  order  to  derive  the  equations  of  motion  of  the  gyrostat,  the  angular  velocities, 
torque,  moments  of  inertia,  and  angular  momentum  of  the  rotor  and  platform  must 
first  be  defined  (5:28). 

•  hi  =  Iiu>i  +  IfU,  =  angular  momentum  oiV  +  H  about  ei 

•  hi  =  liiJi  =  angular  momentum  otV  +  K  about  e,,  i=2,3 

•  ha  =  hi/jJa  +  Wi)  =  angular  momentum  of  Tl  about 

•  /,  =  moment  of  inertia  oiV  •\-'R  about  e„  i  =  1,2,3 

•  I,  =  axial  moment  of  inertia  of  P  about  ei 

% 

•  Ip  =  Ii  —  It  =  axial  moment  of  inertia  of  V  about  ci 

•  Wj  =  angular  velocity  of  V  about  ci,  i=l,2,3 

•  14;,  =  angular  velocity  of  H  about  Ci  relative  to  V 
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Figure  2.1.  Model  of  axial  gyrostat.  The  asymmetric  platform  is  denoted  as  V 
while  the  axisymmetric  rotor  is  denoted  as  Tt.  The  body  frame  Tp  is 
aligned  with  the  principle  axis  and  rotates  relative  to  the  inertial  frame 
Ti.  The  total  angular  momentum  vector  is  denoted  as  h  and  is  constant 
in  direction  and  magnitude. 
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•  ga  =  torque  applied  by  P  on  71  about  cj 

•  c,  =  principal  axes  o(P  +  i  —  1,2,3 

•  f  =  time 

Note  that  the  angular  velocity  vector  is  relative  to  an  inertial  reference  frame  and 
that  the  transverse  moments  of  inertia  of  the  rotor  are  embodied  in  /}  and  I3.  A 
motor  is  incorporated  in  the  gyrostat  for  the  spinup  maneuver  to  provide  torque  to 
both  bodies.  While  the  transverse  torques  between  the  two  bodies  are  present,  they 
are  not  required  for  the  equations  of  motion.  As  noted,  ha  is  the  angular  momentum 
of  the  H  about  ei.  Therefore  the  time  derivative  of  this  angular  momentum  is  the 
axial  torque  applied  by  the  platform  on  the  rotor. 

dha 


The  angular  momentum  components  of  the  gyrostat  may  now  be  written. 


h  =  hj>  + 


(2.3) 


Ax 

/i 

0 

0 

Wi 

I, 

0 

0 

u;. 

h  = 

Aa 

= 

0 

h 

0 

+ 

0 

0 

0 

0 

A3 

0 

0 

h 

W3 

0 

0 

0 

0 

h  = 


/iWi  +  IfU, 
/jWa 

/3U/3 


(2.4) 


(2.5) 


In  order  to  get  the  time  rate  of  change  of  the  angular  momentum  of  the  gyrostat  in 
an  inertial  frame,  the  derivative  of  the  angular  momentum  vector  is  taken  via  the 
standard  convention. 


'd 

+  xh 
dt 


(2.6) 


2-3 


This  yields 


hi  —  waha  +  waha 

0 

ha  +  Wa^i  ~  wiha 

= 

0 

ha  —  Wahl  -f-  wiha 

0 

(2.7) 


The  derivative  of  the  angular  momentum  vector  is  set  equal  to  zero  since  the  gyrostat 
is  modelled  in  an  external  torque  free  environment.  From  Eq.  2.5  we  see  hi  = 
/iwi  +  which  can  be  solved  for  a>i  to  yield 


Wi 


hi  — 

h 


(2.8) 


By  using  the  moment  of  inertia  of  V  about  ei  {Ip  =  /i  —  /,)  and  the  angular 
momentum  of  H  about  e\  {ha  =  /<(a>j  +  u;i))  Ex^.  2.8  can  be  rearranged  to  yield 


CJi 


hj  —  hg 

Ip 


(2.9) 


The  components  of  can  now  be  written  in  terms  of  angular  momentum  and 
moments  of  inertia.  From  Eq.  2.7 


u^h2  —  fjJ^hz 
U^/‘2U^2  —  ^2l^Z 
{I2  ~  /3)c</2<*^ 


finally 


Similarly 


using  Eq.  2.9 


{I2  -  h) 
hh 


h2h3 


— t*^hi  +  ui  ha 


—Wahl  + 


(2.10) 


(2.11) 


(2.12) 
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(2.13) 


and 


dt 


jhz  = 


(2.14) 


- 


h-h 


(2.15) 


hip  *  '  L 

Equations  2.2,  2.10,  2.13  and  2.15  represent  the  complete  system  of  equations  that 
describe  the  spinup  dynamics  of  the  g)rro8tat.  The  cone  angle  (also  known  as  the 
nutation  angle)  is  the  angle  between  the  angular  momentum  vector  and  the  symmetry 
axis  of  the  gyrostat.  In  terms  of  Euler  angles  it  is  defined  as  (5:28) 


V  = 


(2.16) 


S.l  Dimensionless  Equations 

In  order  to  simplify  the  equations  of  motion  further,  the  variables  are  non- 
dimensionalized  using  the  following  transformations  first  developed  by  Guelman 
(4:111)  and  refined  by  Hall  (5:30). 


•  X,  =  hifh  -  dimensionless  component  of  angular  momentum,  *  =  1,2,3 

•  H  =  hal h  -  dimensionless  angular  momentum  of  Ti  about  e\ 

•  t  =  htj Ip  -  dimensionless  time 

•  fif  =  {9alp)lh^  -  dimensionless  torque  applied  by  the  V  to  the  "R,  about  Ci 

The  non-dimensional  equations  will  now  be  with  respect  to  the  non-dimensional 
time  t  and  denoted  as  ()  =  d{)fdt.  Further,  the  moments  of  inertia  are  non- 
dimensionalized  using  the  following  definition  by  Hall  (5:30) 


(2.17) 
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Simple  substitution  into  Eqs.  2.10  2.13  and  2.15  yields  the  following  set  of  dimen¬ 
sionless  equations  on  which  the  rest  of  the  analysis  follows 


Xi 

=  (*2  -  *3)X2X3 

(2.18) 

Xj 

=  (»3Xi  -  p)X3 

(2.19) 

X3 

=  -(*axi-p)x2 

(2.20) 

/i 

=  9 

(2.21) 

Note  that  there  has  been  no  change  in  the  angular  momentum  integral  (Eq.  2.1) 
which  now  becomes 

Xj  +  Xj  +  X3  =  1  (2.22) 

Also  note  that  the  cone  angle  is  reduced  to 

1]  =  cos~^(xi)  (2.23) 


2.2  SpinUp  Boundary  Conditions 


The  spinup  problem  in  this  analysis  usually  begins  with  the  gyrostat  operating 
such  that  the  relative  velocity  between  the  platform  and  rotor  is  close  to  zero.  This 
would  be  the  case  in  both  the  initial  satellite  deployment  or  a  satellite  that  has 
degraded  to  a  flat  spin  for  an  extended  period  of  time.  The  spinup  maneuver  begins 
when  a  small  constant  torque  g  is  applied  to  increase  the  angular  momentum  of  the 
rotor  to  a  value  near  p  =  1.  The  value  of  ^  =  1  is  apparent  if  you  recognize  that  for 
an  inertially  fixed  platform,  u;i  (angular  velocity  of  V  about  ei)  =0.  Substituting 
p  into  Eq.  2.9  yields 


Wi 


Ip/h 


(2.24) 


If  u;i  =  0  then  /z  =  1  for  the  spun  up  condition  (5:58-59). 
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Since  ii  =  g  implies  /x  gt  +  po,  the  lower  boundary  for  n  is  not  necessarily  0. 
Having  Po  /  0  implies  that  the  rotor  is  spinning  relative  to  the  platform  with  some 
initial  angular  velocity  condition.  This  only  affects  the  amount  of  time  required  to 
complete  the  spinup  maneuver  since 


^final  — 


final  —  t^O  1  -  Mo 


(2.25) 


8.3  The  Dynamical  Shape  of  the  Gyrostat 

(5:33-35)  It  can  be  seen  from  Ex).  2.17  that  the  three  dimensionless  inertia 
parameters  retain  their  relative  values  with  respect  to  the  origin2d  moments  of  inertia. 
Mathematically, 

/j  >/*'<=»*,  >»*;  i,  A:  =  1,2,3  (2.26) 

Further,  the  sign  of  and  13  are  determined  by  the  relationship  between  Ip  2uid  /} 
or  I3. 


Ip>Ik  ik<0 

Ip  <  Ik  **  >  0 

ib  =  2,3 

Since  only  ^<1  h  appear  in  £)qs.  2.18  -  2.21,  these  parameters  define  the  dynamical 
shape  of  the  gyrostat.  Using  the  standard  definitions  for  satellite  shape,  the  following 
possibilities  are  available. 

•  if  *3  <  *2  <  0  ^  Oblate  gyrostat 

•  if  *2  >  13  >  0  =►  Prolate  gyrostat 

•  if  i2  >  0  >  is  ^  Intermediate  gyrostat 

The  assumption  is  made  that  12  >  13.  The  other  non-dimensional  inertia  parameter 
simplifies  to  ii  =  I,/Ii  which  is  the  ratio  of  the  axial  inertia  of  H  to  the  axial  inertia 
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Figure  2.2.  Allowable  values  of  and  ta  for  a  gyrostat  with  a  rod  shaped  rotor 
(I'l  =  0).  Oblate  gyrostats  have  negative  values  of  and  is.  Prolate 
gyrostats  have  positive  values  of  is  and  13. 

otV  +  %.  This  implies  that  0  <  I'l  <  1  which  has  the  limiting  physical  shape  of  a 
rod  shaped  rotor  (ti  =  0)  to  a  rod  shaped  platform  (ti  =  1).  The  effect  of  t’l  on  the 
gyrostats  dynamics  are  described  in  Hall’s  work  (5:33-35).  He  shows  that  by  using 
the  triangle  inequalities  for  the  moments  of  inertia,  limits  the  physically  possible 
values  of  is  and  13  through  the  following  relations. 

(1  -  2ix  +  iiis)l{i2  ~  m)  (/) 

*3  =  ■  (213  -  ills  -  1)1  {is  -  *i)  {II) 

{l-iii2)f{2-ix-i2)  {III) 

A  graphical  presentation  is  shown  in  Figure  2.2  for  a  gyrostat  that  has  ii  =  0.  In 
accordance  with  the  relations  given  above,  oblate  gyrostats  will  have  negative  is  2und 
is  values  while  prolate  gyrostats  will  have  positive  values  of  is  and  is.  Intermediate 
gyrostats  will  not  be  discussed  in  this  thesis. 


Hall  also  demonstrates  that  ti  impacts  the  dynamics  of  the  gyrostat  through 
its  effect  on  the  initial  conditions.  Note  that  the  definition  of  oblate  and  prolate  do 
not  involve  t'l.  Because  spinup  typically  begins  with  the  platform  and  rotor  spinning 
as  a  single  rigid  body  (a>«  =  0),  the  following  relations  hold. 


hi  —  IiLJ\  + 

ha  =  IaU>t  +  /jU>i  =  laUl 


Since 


and 


(2.27) 


(2.28) 


using  Eq.  2.27 

fji  = 

ft  ii  ii 

This  initial  spin  can  be  about  either  the  major  or  minor  moment  of  inertia  axis  of 
the  gyrostat.  For  example,  an  aJl  spun  prolate  spacecraft  can  be  represented  by  two 
possibilities  (6:644). 


•  if  ii  >  12  >  *3  >  0  fben  Ci  is  the  major  axis  and  the  state  is  Xi  =  1,  X2  =  X3  =  0, 

/i  =  l. 

•  if  12  >  *1  >  *3  >  0  then  62  is  the  major  axis  and  the  state  is  X2  =  1, 
Xi  =  X3  =  /i  =  0.  This  is  commonly  referred  to  as  a  fiat  spin  condition. 

The  reader  is  referred  to  reference  (5)  for  an  in-depth  discussion  of  the  dimensionless 
parameters. 
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8.4  Spinup  Problems 

To  describe  the  nature  of  the  different  spinup  problems,  an  understanding  of 
the  system  dynamics  when  =  0  is  first  required.  As  can  be  seen  from  Eq.  2.21, 
setting  9  =  0  implies  that  the  angular  momentum  of  H  about  e\  is  a  constant  value. 
Because  the  angular  momentum  first  integral  cannot  be  violated  in  a  torque  free 
environment,  the  motion  of  the  gyrostat  is  confined  to  one  dimensional  curves  on 
the  momentum  sphere  defined  by  Eq.  2.22.  Hall  (5:Chapter  5)  demonstrated  that 
when  9  =  0,  Eqs.  2.18-2.20  reduce  to  a  single  integral  which  is  solved  using  Jacobi’s 
elliptic  functions.  In  Hall’s  development  of  the  reduction  of  quadrature  for  9  =  0, 
he  derives  an  energy  like  constant,  y,  which  defines  lines  of  constant  “energy”  on 
the  momentum  sphere  (5:39-42).  The  “energy”  term,  y,  represents  a  functional 
combination  of  the  three  integrals  of  motion  when  9  =  0.  The  three  integrals  are 
total  angular  momentum,  axial  angular  momentum,  and  rotational  kinetic  energy. 

y-\  {(*2  +  *3)^1  “  (*2  “  *3)(^2  “  2:3)  -  Apxi  -  (12  -h  13)}  (2.30) 

Because  this  constant  is  only  valid  for  ft  =  constant,  the  momentum  sphere  has  a 
different  topology  for  each  value  of  p.  Figure  2.3  depicts  a  progression  of  momentum 
spheres  for  a  near  axisymmetric  gyrostat  as  p  increases.  This  is  essentially  a  graphical 
presentation  of  a  spinup  maneuver  for  either  an  oblate  or  prolate  gyrostat.  In  Section 
2.3,  we  discussed  the  non-dimensional  moments  of  inertia  Z2  and  Z3.  If  both  these 
values  are  negative,  the  gyrostat  is  oblate  and  its  energy  for  a  particular  value  of  p 
is  associated  with  a  polhode  that  lies  in  the  upper  (northern)  half  of  the  momentum 
sphere  (note:  “upper”  and  “lower”  are  in  reference  to  the  spheres  shown  in  Figure 
2.3  regardless  of  the  coordinate  axis).  If  12  and  Z3  are  both  positive,  the  gyrostat  is 
prolate  and  its  energy  is  associated  with  polhodes  that  lie  in  the  lower  (southern) 
half  of  the  momentum  sphere.  Note  that  the  number  of  equilibrium  points  on  the 
momentum  sphere  vary  as  the  values  of  p  change.  Hall’s  work  demonstrates  that 
for  p  less  than  the  smaller  of  |z2|  or  (I’sl,  there  are  six  equilibrium  points.  The  six 
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Table  2.1.  Equilibrium  Points 


Equilibrium  Pis 

Validity 

(±1,0,0) 

Vg 

1 

g=0 

1 

g=o 

equilibrium  points  are  listed  in  Table  2.1  and  were  first  published  in  Guelman’s  work 
(4:112). 

The  points  and  P^^  in  Figure  2.3  are  centers  for  an  oblate  and  prolate 
gyrostat  respectively.  The  two  points  indicate  the  unstable  equilibria  for  unstable 
flat  spin  motion.  The  two  points  labeled  located  within  the  unstable  flat  spin 
separatrix  correspond  to  the  centers  of  stable  flat  spin  motion.  As  /i  increases, 
the  separatrices  associated  with  flat  spin  motion  migrate  to  the  prolate  equilibrium 
point.  At  n  equal  to  the  smaller  of  |i2|  or  |t3|,  the  two  saddles  of  the  unstable  flat 
spin  motion  converge  to  form  a  single  saddle  at  the  south  pole  and,  as  a  result, 
with  the  smaller  of  |i2|  or  |i3|  <  n  <  the  larger  of  1*21  or  |*3|,  there  are  only  four 
equilibrium  points.  Finally,  at  p  equal  to  the  larger  of  |*2|  or  |*3|,  the  two  stable 
flat  spin  centers  converge  with  the  saddle  at  the  south  pole  to  form  a  center.  For 
fi  greater  than  the  larger  of  |*2|  or  |*3|  only  two  equilibrium  points  remain  and  no 
further  bifurcations  occur.  It  is  important  not  to  visualize  the  polhodes  as  circular 
disks.  The  momentum  sphere  contains  separatrices  which  affect  all  the  polhodes 
on  the  sphere.  A  more  realistic  concept  of  what  is  seen  as  circles  is  a  surfeice  of  a 
hyperboloid  or  a  “Pringles  potato  chip”. 

Unfortunately,  the  spinup  problem  involves  the  case  when  g  ^0.  When  5  #  0, 
there  are  no  polhodes  defined  by  Eq.  2.30  since  pi  is  no  longer  constant.  Further,  the 
separatrices  which  separate  the  different  kinds  of  motion  in  a  y  =  0  system  no  longer 
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separate  different  regions  of  phase  space.  It  is  entirely  possible  for  a  trajectory  to 
cross  the  “instantaneous”  separatrix  as  /i  is  increasing.  The  term  “instantaneous” 
is  used  since  at  each  value  of  n,  there  exists  a  range  of  energy  values  which  can  be 
plotted  as  instantaneous  polhodes  on  the  momentum  sphere.  One  way  to  envision 
the  dynamics  is  to  think  of  the  trajectory  moving  on  the  momentum  sphere  as 
the  instantaneous  separatrices  migrate  to  P^.  An  important  conceptual  point  is  to 
realize  that  for  near  axisymmetric  gyrostats,  the  trajectory  will  not  move  towards 
the  oblate  or  prolate  equilibrium  points.  In  fact,  without  a  significant  asymmetry 
in  the  transverse  moments  of  inertia,  the  xi  component  of  angular  momentum  will 
stay  at  approximately  the  same  value.  This  implies  that  the  cone  angle  at  the  end 
of  spinup  is  not  going  to  vary  significantly  from  .  initial  value. 

As  mentioned  earlier,  trajectories  for  oblate  gyrostats  will  begin  in  the  northern 
hemisphere.  The  instantaneous  separatrices  on  the  momentum  sphere  will  migrate 
toward  the  south  pole,  so  the  oblate  trajectory  will  not  experience  any  separatrix 
crossings.  Prolate  gyrostats  begin  spinup  in  the  southern  hemisphere.  As  the  spinup 
progresses,  the  trajectory  is  influenced  by  the  instantaneous  separatrices.  Specifi¬ 
cally,  as  the  trajectory  crosses  an  instantaneous  separatrix,  its  frequency  decreases 
until  the  second  bifurcation  occurs.  At  this  point,  the  trajectory  paths  look  similar 
to  those  present  before  the  separatrix  crossing  and  the  frequency  should  increase 
again.  Transverse  spinup  of  a  prolate  gyrostat  begins  spinup  near  one  of  the  stable 
flat  spin  equilibrium  points.  In  order  to  keep  the  trajectory  near  the  equilibrium 
point,  the  gyrostat  must  have  an  extremely  small  torque  or  substantial  asymme¬ 
try.  Since  this  thesis  deals  with  near  axisymmetric  gyrostats,  the  second  option  is 
not  considered.  If  the  gyrostat  has  too  much  torque,  the  trajectory  will  cross  the 
instantaneous  c  iparatrix.  Obviously,  with  a  higher  torque,  this  will  happen  quicker. 

Motion  near  a  separatrix  is  highly  non-linear  and,  therefore,  difficult  to  approx¬ 
imate.  As  a  result,  we  expect  to  experience  some  difficulties  associated  with  prolate 
and  transverse  spinup.  As  mentioned  above,  the  problem  of  defining  exactly  where 
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the  instantaneous  separatrix  crossing  occurs  is  extremely  difficult.  The  reason  is  the 
constantly  changing  nature  of  the  instantaneous  separatrices  as  the  spinup  maneuver 
progresses.  To  define  the  separatrix  crossing,  the  rate  at  which  the  instantaneous 
separatrices  are  moving  on  the  momentum  sphere  along  with  the  exact  position  and 
rate  of  movement  of  the  angular  momentum  vector  on  the  momentum  sphere  would 
have  to  be  known.  Given  the  extremely  non-linear  nature  of  the  motion  in  this  region 
along  with  an  infinite  number  of  initial  conditions;  the  difficulty  of  the  problem  can 
be  appreciated. 

2.5  Summary 

In  this  chapter,  we  have  derived  the  non-dimensional  equations  of  motion  that 
describe  the  spinup  of  an  axial  gyrostat.  The  spinup  boundary  conditions  have  been 
specified  and  a  discussion  of  the  three  types  of  spinup  problems  in  terms  of  their 
trajectories’  relationship  to  the  nonlinear  effects  of  separatrices  has  been  presented. 
In  the  next  chapter,  the  straightforward  expansion  method  is  used  to  outline  the 
general  flow  of  an  approximate  solution  for  the  governing  equations. 
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III.  The  Straightforward  Expansion 

We  begin  the  analysis  using  the  straightforward  expansion  method  as  a  pos¬ 
sible  perturbation  technique.  The  rationale  behind  using  this  approach  is  twofold. 
First,  by  employing  the  straightforward  expansion  method,  a  general  flow  of  the 
solution  using  perturbation  theory  is  presented.  Second,  the  method  is  employed  to 
determine  if  valid  solutions  can  be  obtained  for  larger  torque  values.  A  well  known 
limitation  of  the  straightforward  expansion  method  is  that  it  produces  terms  which 
are  proportional  to  the  independent  variable,  t.  These  terms,  commonly  referred  to 
as  secular  terms,  invalidate  an  approximate  solution  at  large  values  of  the  indepen¬ 
dent  variable  because  they  grow  without  bound.  However,  for  larger  values  of  g, 
Eq.  2.25  shows  that  the  amount  of  time  required  for  the  spinup  problem  is  reduced. 

We  wish  to  derive  an  approximate  solution  for  the  spinup  of  a  near  axisym- 
metric  gyrostat  which  have  the  following  equations  of  motion  (Eqs.  2.18-2.21). 


Xi  =  (t3-t3)X2X3 

Xj  =  (»3Xi  -  n)X3 

X3  =  -(tjXi  -  li)X2 

it  =  9 


This  system  is  subject  to  the  following  initial  conditions. 


The  small  parameter  which  we  use  in  the  approximation  is  the  difference  in  transverse 
inertias 

«2  -  »3  =  ae  (3.2) 

where  a  is  (?(1).  We  begin  with  the  assumption  that  each  component  of  amgular 
momentum  has  a  solution  of  the  form 


Xi  =  Xio  +  £Xii  +  £^Xi2  +  O(c^);  1  =  1, 2, 3  (3.3) 

Further,  the  time  derivative  of  each  component  takes  the  form 

Xi  =  Xio  +  £Xii  +  £^Xi2  +  O(e^);  i  =  1, 2, 3  (3.4) 

Note  that  Xj  refers  to  the  total  component  of  angular  momentum  while  Xij  refers 
to  the  terms  which  make  up  the  solution  for  x^.  Substituting  Eqs.  3.3  and  3.4  into 
Eqs.  2.18-2.20  and  collecting  terms  of  e  yields  the  following  system  of  equations. 


0(1): 

iio  = 

0 

(3.5) 

X20  = 

-(M  -  *3Xlo)*30 

(3.6) 

^30  = 

in  -  *3*lo)*20 

(3.7) 

0{£): 

ill  = 

a  X20*30 

(3.8) 

iji  = 

-(/I  -  i3Xlo)*31  +  *3*11*30 

(3.9) 

*31  = 

(/*  —  i3*lo)*31  ~  0*10*20  —  *3*11*20 

(3.10) 

Notice  the  important  contribution  of  selecting  the  difference  in  transverse  in¬ 
ertias  as  the  small  p2U'ameter.  Elquations  2.18-2.20  are  transformed  from  a  set  of 
non-linear  differential  equations  to  systems  of  linear  differential  equations  in  each 
power  of  £.  Solving  the  0(1)  system  of  equations  first,  we  immediately  see  that 
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Xio  =  constant,  leaving  a  set  of  coupled  linear  differenti2J  equations  of  the  form 


X  =  yl(t)X 


(3.11) 


where 


A(t)  = 


0  -(/<  -  13X10) 

(/i  -  *3X10)  0 


ti  =  9t  +  IMi 


(3.12) 

(3.13) 


This  system  has  the  solution 

X  =  *(t,0)Xo  (3.14) 

where  Xo  are  the  initial  conditions  and  $  is  the  state  transition  matrix.  While 
coupled  systems  of  differential  equations  with  time  varying  coefficients  are  typically 
difficult  to  solve,  this  system  of  equations  contiuns  a  simple  property  which  permits 
a  simple  solution.  We  make  use  of  the  following  relation  (7:p.  600) 


If  A(t\)A{t2)  —  A(<2)i4(<i);  then  ♦(t,T)  =  exp  j  A{<r)da  (3.15) 


The  A  matrix  in  this  system  of  equations  satisfies  the  commutative  property  and  so 
the  integration  is  easily  done: 


+  Po  —  t3Xio)do  =  ^  +  (po  —  isXiojt 


(3.16) 


We  now  make  the  following  definitions. 


k  =  Po  —  *33^10 

(3.17) 

m  =  Y + “ 

(3.18) 
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The  resulting  state  transition  matrix  is 


0)  =  exp 


0 

-m ' 

./(<) 

0 

cos  /(<) 
sin  f{t) 


—  sin  f(t) 
cos  f{t) 


(3.19) 


Note  that  if  we  were  to  stop  the  expansion  right  here,  we  would  be  deriving  an 
exact  solution  to  the  equations  of  motion  for  a  completely  axisymmetric  g}rrostat 
(i.e.  e  =  0).  Application  of  the  initial  conditions  to  the  0(1)  system  solutions  yields 
the  following  result. 


Xio  =  Xi(0) 

(3.20) 

X20  =  X2(0)cos/(t)  —  X3(0)sin/(t) 

(3.21) 

X30  =  X2(0)  sin  /(<)  +  X3(0)  cos  f{t) 

(3.22) 

We  now  proceed  to  the  0(ff)  system  of  equations  shown  here  again  for  conve¬ 
nience. 


ill  =  a  X20X30 

iji  =  — (^  —  *3aJio)x3i  -f-  133:11X30 

X3I  =  (a*  —  *3Xio)X2i  —  aXioX20  —  *3X11X20 


Starting  with  the  xn  term 


Xll  = 


/  {0x20X30}  dar  (3.23) 

Jo 

jf  a  {x2(0)x3(0)  cos^  /(cr)  -|-  (x|(0)  -  X2(0))  cos  f{a)  sin  /(<t) 
-X2(0)x3(0)  sin^  /(<^)}  da  (3.24) 
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The  result,  using  the  symbolic  manipulation  program  Matkematica  2.0  for  SPARC 
by  Wolfram  Research  Inc.  (15),  is 


xii  =  S 


-  C 


-I-  gt) 


v/? 


6i+C 


y/9 


bi-S 


*2 


with 


(3.25) 


5i  =  a,f^(x2i0)^cos- - X3(0)^  cos  —  +  2x2(0)a:3(0)8in  — 

Y  0^  V  9  9  9 ) 

62  =  oi\ f — iC2(0)^  sin  —  +  X3(0)^  sin  —  +  2x2(0)x3(0)  cos  — ^ 

9  9  9  } 


The  Fresnel  Sine,  denoted  by  “S”,  and  Fresnel  Cosine,  denoted  by  “C”,  are  functions 
that  are  defined  as  (1:300) 


C(x)  = 
5(x)  = 


dt 

dt 


Figure  3.1  depicts  both  of  these  functions.  Note  that  as  t  — »  00,  both  functions 
cisymptotically  approach  5.  The  Fresnel  Sine  and  Fresnel  Cosine  can  be  approxi¬ 
mated  using  the  following  functions. 


C(x)  = 
S(x)  = 


i  +  psin(p) 


where 


1  -I-  .926x 

2  +  1.729X  -I-  3.104x2 


(3.26) 
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Figure  3.1.  Fresnel  Sine  and  Fresnel  Cosine  Functions. 


1 


2  +  4.124*  +  3.429*2  +  6.670*3 


(3.27) 


These  approximations  are  valid  to  2  x  10“®  for  0  <  *  <  oo  (1:302).  An  important 
consideration  is  the  singularity  at  *  =  —0.4932  in  Eq.  3.27.  This  inconvenience  can 
be  accommodated  by  using  the  following  symmetry  relations  (1:301). 

C(-*)  =  -C{x) 

5(— *)  =  — 5(*) 

With  *11  determined,  we  now  proceed  to  the  set  of  coupled  differential  equations  for 
*21  and  *31.  Eqs.  3.9  and  3.10  form  a  system  that  has  the  form 


X  =  A{t)X+Bu{t) 


(3.28) 
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where  A(t)  is  the  same  as  Eq.  3.12,  B  is  the  identity  matrix  and 


^30*3^^11 

—  QXxqXjq  —  *3X11X20 

This  is  simply  the  C?(l)  system  of  equations  (cf.  Eq.  3.11)  with  a  forcing  function 
u.  This  system  has  a  solution  of  the  form 


u(0  = 


X(t)  =  <&(<,0)X(0)  +  $(<,0)  /‘  $-*(a)u(a)d<T  (3.30) 

Jo 

The  state  transition  matrix  $  is  defined  by  Eq.  3.19.  The  first  term  of  the  solution 
must  be  set  to  zero  since  the  initial  conditions  were  matched  in  the  0(1)  system  of 
equations.  The  multiplication  of  the  integrand  results  in 


3:3*3X11  -  I  sin(2/)a:ioa  X2  +  5x3a  Xio  -  5  cos(2/)x3Q  xio 
— X2*3Xii  +  5  sin  (2/)xioa  X3  -  ^X2a  xio  ~  5  cos  (2/)x2a  Xio 


(3.31) 


Recall  that  /  =  /(f)  which  is  defined  by  Eq.  3.18  and  Xi  are  the  initial  conditions. 
The  first  and  third  terms  in  each  row  of  the  integrand  will  produce  secular  terms, 
but  the  second  and  fourth  terms  will  not.  Integrating  each  row  of  Eq.  3.31  and 
multiplying  the  resulting  matrix  by  $  yields  the  following  results  (secular  terms  are 
over-bracketed).  The  terms  ai,  and  5*;  k  =  1,2, 3,4  are  located  in  Appendix  A. 
Further,  the  terms  X3  and  X3  are  defined  as  X2(0)  and  X3(0)  respectively. 


X21  =  xiofx3  cos  (/)  -^a  izXz{x\  -  x|)  cos  (/) 
I  ig 


-I-  *3X3(x2  -  X3)  cos  (/)  cos  2/  +  xiofx3  sin  (/) 

L 

-  —a  *3X2(x^  -  x\)  sin  (/)  -  *3X2(x|  -  x^)  cos  (2/)  sin  (/) 

4flf  4jf 


5 

k\j2lir 

■al-k-C 

kyflfi( 

•a2 
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(3.32) 


+  5 


y/^{k  +  gt) 

y/^{k  +  gt) 

V?  J 

•  uO  T  ^ 

[  ^ 

a4 


-  ttQ!  *3X2X3  cos  if)  sin  (2/)  -  13X3X3  sin  (/)  sin  (2/) 

zg  Ig 


X31  =  -  ixioxj  cos  (/) -^a  i3X2(x|  -  x^)  cos  (/) 
2  4<7 


+  *3X2(x3  -  Xj)  COS  (/)  COS  (2/)  +  fxioX3  sin  (/) 

45  2 

-  —a  t3X3(x2  -  x|)  sin  (/)  +  —  a  *3X3(x2  -  Xg)  cos  (2/)  sin  (/) 
4flf  4flf 


s 

kyj^ 

•  61  +  C 

L  ^  \ 

[  V?  _ 

+  S 


yj^{k  +  gt) 


y/9 


63  +  C 


62 

+  gt) 


V9 


•64 


+  :;^oc  *3X2X3  cos  (/)  sin  (2/)  -  *3X2X3  sin  (/)  sin  (2/)  (3.33) 

^5  2<? 


Notice  that  there  are  no  secular  terms  in  the  xn  solution  while  they  do  appear  in 
the  X21  and  X31  solutions. 


3.1  Straightforward  Expansion  Results 

Because  we  know  that  secular  terms  are  present  in  the  solution,  we  test  the 
solution  on  an  oblate  gyrostat  to  avoid  any  problems  with  instantaneous  separa- 
trix  crossings.  Figures  3.2  through  3.5  depict  the  validity  of  this  solution  with 
the  numerical  integration  of  the  equations  of  motion.  The  numerical  integration 
is  done  using  a  4**  order  Runge-Kutta  subroutine  in  the  computer  program  Mat- 
lab  4-Oa  by  the  Mathworks  Inc.  (10).  The  initial  conditions  will  always  be  written 
at  the  top  of  the  graph  as  the  “state  vector.”  The  state  vector  is  interpreted  as 
[xi(0);x2(0);x3(0);/t(0)].  Further,  the  x  axis  is  represented  as  p  and  is  essentially  a 
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time  axis  from  the  beginning  to  the  end  of  spinup.  As  can  be  seen,  the  solution  for 
the  Xi  component  of  angular  momentum  is  very  good.  However,  as  the  amount  of 
time  required  for  spinup  increases  with  smaller  values  of  g,  the  approximate  solution 
begins  to  get  out  of  phase  with  the  numerical  solution  (see  Figure  3.3).  The  X2 
component  of  angular  momentum  captures  the  increasing  frequency  of  the  spinup 
maneuver  but,  2is  g  decreases  the  secular  terms  become  dominant  and  invalidate  the 
solution.  Figure  3.6  and  3.7  depict  the  maximum  error  seen  throughout  spinup  for 
different  t final  or  g  values.  The  xj  and  X3  components  of  angular  momentum  show  a 
large  divergence  in  error  as  the  spinup  times  increases.  As  a  result,  this  approximate 
solution  is  invalid  for  ^  <  0.1.  The  xi  component  of  angular  momentum  displays 
very  little  error  as  the  spinup  time  is  increased.  This  is  due  to  the  fact  that  there  Me 
no  secular  terms  in  the  xi  solution.  It  is  interesting  to  note  that  the  Xi  approximate 
solution  does  capture  the  decreasing  amplitude  as  spinup  time  increases.  The  error 
seen  in  Figure  3.6  is  due  to  the  out  of  phase  condition  that  develops.  However,  be¬ 
cause  the  scale  of  the  Xi  component  is  so  small,  (i.e.  increments  as  small  as  0.0005), 
the  error  would  still  be  less  then  one  percent  even  if  the  approximate  solution  were 
completely  out  of  phase  with  the  numerical  solution.  For  smaller  vadues  of  g  (for 
example  g  =  0.001),  the  approximate  solution  does  get  completely  out  of  phase 
with  the  numerical  solution  about  half  way  through  the  spinup  maneuver.  However, 
when  fi  —*  I,  the  approximate  solution  is  back  in  phaise  with  the  numerical  solution. 
Further,  the  amplitude  error  is  minimal  throughout  spinup.  With  such  small  error 
in  the  Xi  component  of  angular  momentum,  the  cone  angle  described  by  Eq.  2.23 
is  useful  at  any  time  during  the  spinup  maneuver.  Note  that  the  example  in  this 
anailysis  began  spinup  relatively  close  to  the  oblate  equilibrium  point  located  at  the 
north  pole  of  the  momentum  sphere.  The  region  of  validity  of  the  Xi  approximate 
solution  will  be  explored  further  in  Chapter  4. 
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Figure  3.2.  Straightforward  expansion  method  approximate  solution  of  the  Xi  com¬ 
ponent  of  angular  momentum  of  an  oblate  gyrostat  with  initial  condi¬ 
tions  xi  =  0.9, xj  =  X3  =  0.30822  and  e  =  0.01;  5  =  0.1  . 


MBW  v«clor<4.9^0S22;30a22;0] 


mu 


Figure  3.3.  Straightforward  expansion  method  approximate  solution  of  the  Xi  com¬ 
ponent  of  angular  momentum  of  an  oblate  gyrostat  with  initial  condi¬ 
tions  xi  =  0.9, X2  =  X3  =  0.30822  and  e  =  0.01;  g  =  0.01. 
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•tu*  wclor  t.9:.30622;.30822:01 


Figure  3.4.  Straightforward  expansion  method  approximate  solution  of  the  Xj  com¬ 
ponent  of  angular  momentum  of  an  oblate  gyrostat  with  initial  condi¬ 
tions  xi  =  0.9, xa  =  X3  =  0.30822  and  e  =  0.01;  5  =  0.1. 


v«ctor>(.0;30e22;3Qe22;0] 


Figure  3.5.  Straightforward  expansion  method  approximate  solution  of  the  xa  com¬ 
ponent  of  angul2ir  momentum  of  an  oblate  gyrostat  with  initial  condi¬ 
tions  xi  =  0.9, xa  =  X3  =  0.30822  and  e  =  0.01;  g  =  0.01. 


3-11 


StM*  Vactof  (.7;.4;.59161  ;0| 


Figure  3.6.  Maximum  error  detected  in  the  Xi  component  of  angular  momentum  as 
a  function  oft  final  (or  g)  for  an  oblate  gyrostat  using  the  straightforward 
expansion  method. 


stm  Vador  |.7;.4;.59iei  ;0) 


Figure  3.7.  Maximum  error  detected  in  the  X2,  X3  component  of  angular  momentum 
as  a  function  of  t final  (or  g)  for  an  oblate  gyrostat  using  the  straight¬ 
forward  expansion  method. 
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3.8  Conclusions 


The  analysis  thus  far  has  demonstrated  that  the  straightforward  expansion 
method  provides  reasonable  results  for  the  small  class  of  oblate  gyrostats  that  have 
a  near  axisymmetric  configuration  with  a  very  large  torque  value.  However,  the 
secular  terms  and  the  length  of  the  equations  in  the  solution  both  contribute  to 
its  general  failure.  The  approach  did  point  out  the  areas  with  which  the  solution 
needs  to  address  itself.  The  equations  of  motion,  when  numerically  integrated,  will 
produce  solutions  which  have  time  varying  frequency  and  amplitude.  The  method 
of  multiple  scales  is  ideally  suited  to  handle  these  problems  and  is  the  subject  of  the 
next  chapter. 
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IV.  The  Method  of  Multiple  Scales 

As  seen  in  Chapter  3  and  specifically  Figure  3.7,  the  straightforward  expansion 
method  only  produces  acceptable  results  for  times  (or  equivalently  torque  values) 
that  are  very  small  (large)  compared  to  1/c.  This  is  an  inherent  problem  with  series 
expansions  or  perturbation  solutions  of  nonlinear  differential  equations.  The  problem 
arises  due  to  the  fact  that  0(e)  and  higher  terms  produced  by  the  straightforward 
expansion  method  are  proportional  to  the  independent  variable,  t.  These  terms, 
commonly  referred  to  as  secular  terms,  destroy  the  fundamental  assumption  of  a 
perturbation  solution.  Specifically,  in  a  series  expansion,  each  successive  power  of 
the  small  parameter,  e,  has  decreasing  importance  in  the  series  solution.  By  being 
proportional  to  the  independent  variable,  these  terms  rapidly  become  large  in  relation 
to  the  (?(1)  term  (11:24-25).  This  result  is  also  known  as  a  nonuniform  expansion. 

The  Method  of  Multiple  Scales  was  developed  to  overcome  this  deficiency. 
Eq.  3.32  contuns  the  0(e)  term  for  the  Xj  component  of  angular  momentum  in  the 
straightforward  expansion  solution.  The  functional  dependence  of  X3  on  t  and  e  can 
also  be  written  as 

X2(t;e)  =  X2(t,et,e*t,---;e)  (4.1) 

Or 

X2(t;  e)  =  x‘2(to,  ti,  t2,  •  •  • ;  e)  (4.2) 

where  the  tn  are  defined  as 


<0  =  h  =  et,  <2  =  (4.3) 

Note  that  the  are  different  time  scales.  If  e  =  ^,  then  variations  on  to  could  be 
observed  on  the  second  hand  of  a  watch.  Variations  on  ti  could  be  observed  on  the 
minute  hand.  Further,  variations  on  <2  could  be  observed  on  the  hour  hand  and  so 
on  (12:122-123).  As  a  result,  the  functional  dependence  of  X2  on  a  single  independent 
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variable,  t,  changes  to  dependence  on  To  incorporate  this  dependence 

in  the  solution  we  use  the  differential  operator 


±_d_ 

dt~  dto^^  dti  at^'^ 
The  assumed  series  solution  for  Xi  is  now 


(4.4) 


X,  =  +  +  *  =  1.2,3  (4.5) 


If  these  relations  are  substituted  into  a  differential  equation  and  powers  of  e  collected, 
they  result  in  a  series  of  partial  differential  equations  as  functions  of  the  time  scales 
tn.  Solutions  to  these  PDE’s  involve  the  determination  of  functions  which  cause 
secular  terms  to  vanish  (2:120-123). 

We  begin  the  multiple  scales  solution  approach  with  the  four  first  order  differ¬ 
ential  equations  that  describe  the  system. 


Xi  =  (t3-i3)X2X3  (4.6) 

X2  =  (t3Xi-/i)X3  (4.7) 

X3  =  -(l2Xi-^)X2  (4.8) 

A  =  5  (4.9) 


This  system  is  subject  to  the  following  initial  conditions 


x(0)  = 


xiW 

X3(0) 

X3(0) 

MO) 


(4.10) 


Note  that  Xi  refers  to  the  total  component  of  angular  momentum  while  Xij  refers  to 
the  terms  which  mjdce  up  the  solution  for  x,.  The  small  parameter  about  which  we 
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want  to  perturb  the  system  is  the  difference  in  the  transverse  inertias 


*2  —  t3  =  a  e  (4.11) 

where  a  is  C7(l).  We  only  use  two  time  scales  and  seek  bounded  solutions  of  the 
form 


=  Xio{to,ti)  +  £Xii{to,ti)  +  £^Xi2{toyti)  +  0(£^);  i  =  1,2,3  (4.12) 


The  fast  time  “to”  is  defined  as  t  while  the  slow  time  “ti”  is  defined  as  e  t.  In  order 
to  incorporate  the  different  time  scales,  the  following  differential  operator  is  used. 


£_d_  _d_ 

dt  dto  ^  ^  dti 


(4.13) 


In  order  to  insure  that  angular  momentum  is  conserved,  we  impose  the  constraint 


l|x||’  =  l,Vt>0 


(4.14) 


Expanding  this  relation  using  E)q.  4.12  we  see 
1  =  ||x(to,  ti;  e)||^  =  X  •  X  =  ||xo||^  +  2eXo  •  Xi  +  e*  [||xi|l*  +  2xo  •  X2]  +  •  •  •  (4.15) 


This  imposes  an  orthogonality  condition  on  the  C>(e)  terms.  Plugging  Eqs.  4.11,  4.12 
and  4.13  into  Eqs.  4.6  through  4.8  and  collecting  terms  of  C7(c)  yields  the  following 
systems  of  equations. 


0(1): 


dxio 

^0 

9X20 

9to 

9X30 

9to 


0 


(*3X10  -  fi)X30 


—(*33710  —  P)3720 


(4.16) 

(4.17) 

(4.18) 
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(4.19) 


Subject  to 


;Eio(0,0)  =  Xi(0),  i  =  1,2,3 
llxof  =  1,V  (>  0 


0(£); 


Oto 

9x21 

dto 

dxzi 

dto 


Subject  to 


dxio 

a  X2oXao  — 


(*3(®10  —  li)X3l  +  *3iCnaC30 


dx^ 

~dU 


—  (t3iClO  —  /*)*21  —  *33Jlia?20  —  O'  ^loJ^ao  — 


dx: 


30 


dti 


(4.20) 

(4.21) 

(4.22) 


x.i(0,0)  =  0,  t  =  1,2,3 
Xo  •  Xi  =  0,  V  t  >  0 


Starting  with  the  0(1)  system  of  equations,  the  xio  term  has  the  solution 


xio(to»<i)  =  •Ai(<i)  (4.23) 

where  i4i(<i)  is  the  “constant”  of  integration.  The  term  “constant”  is  used  to  indicate 
that  although  the  term  has  time  dependance  in  ti,  it  is  constant  with  respect  to  to. 
In  order  to  determine  the  dependence  of  ti  on  zio,  we  must  solve  the  next  higher 
order  of  the  system  of  equations,  0(e).  The  xjo  and  X30  terms  are  coupled  partial 
differential  equations.  Note  that  p  is  only  a  function  of  the  fast  time  to.  This  can  be 
seen  by  noting  that  fi  reflects  the  energy  that  is  being  put  into  the  system.  It  is  the 
result  of  a  simple  ordinary  differential  equation  and,  as  a  result,  has  no  dependence 
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on  the  slow  time.  Making  use  of  E)q.  4.23,  the  system  of  equations  can  be  written  as 


=  (i3Ai(ti)  -  n{to)]  X3o{to,  ti) 

-MM]a;2o(to,ti) 


(4.24) 

(4.25) 


This  system  can  be  represented  as 


— Xo  =  B{to,ti)XQ{tQ,ti) 

UlQ 


X>  = 


X2j 

^3i 


and  the  B  matrix  is 


B 


0  [t3A(ti)  -  /i(to)] 

-  [i3A(ti)  -  tl{to)]  0 


(4.26) 

(4.27) 


(4.28) 


This  matrix  is  very  similar  to  the  matrix  developed  in  Chapter  3  using  the  straight¬ 
forward  expansion  method.  It  still  satisfies  the  property  of  commutativity  and  so  we 
make  use  of  the  following  relation  (7:600) 

If  B{ti)B{t2)  =  B{t2)B{ti)]  then  $(t,r)  =  exp  J  B{a)d<T  (4.29) 
and  the  system  of  equations  has  the  solution 


X  =  $(f,0)Xo(0,ti) 


(4.30) 


The  matrix  $(t,  0)  is  termed  the  state  transition  matrix  and  Xo  are  the  constants 
(with  respect  to  to)  of  integration  of  the  system.  The  integration  of  the  B  matrix  is 
easily  done 

/  {gtr  +  /io  -  izA{ti))  d<r  =  -^  -I-  (/io  -  t3A(<i))  <o  (4.31) 

Jo  e, 
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To  simplify  the  equations  we  make  the  following  definitions. 


fi^o)  =  ^  +  kto 


As  a  result,  the  state  transition  matrix  is  again 
$  (t,  0)  =  exp 

Finally,  the  solution  to  the  set  of  coupled  PDF’s  (Eqs.  4.17  2uid  4.18)  is 


0 

-fito)  ' 

cos /(to)  -sin /(to) 

.  /(to) 

0 

sin /(to)  cos /(to) 

(4.32) 

(4.33) 


(4.34) 


X7oiU,ti)  =  A‘i(ti)cosf{to)-A3{ti)sinf{to)  (4.35) 

X3o(<o,«i)  =  A2(<i)sin/(to)  +  A3(ti)cos/(to)  (4.36) 


We  now  proceed  to  solve  the  0(e)  system  of  equations  to  find  the  <i  dependence 
in  XiQ.  The  (!7(e)  equations  are  listed  here  again  for  convenience. 


_ 

^0 

dX2l  _ 
dtf) 

dx3i  _ 

9to 

Noting  Eq.  4.23,  the  Xn 


dx 


a  X20X30  — 


10 


dti 


(taario  —  ti)x3i  +  *3X11X30  — 


dx 


20 


dti 


—  (*3®10  ~  ~  *3*113^20  —  01  X1QX20  — 


term  can  be  written 


dx3o 

dti 


(4.37) 

(4.38) 

(4.39) 


^  =  -A;(ti)  +  a  X 

[A2(ti)  cos  /(to)  -  A3(<i)  sin  /(to)]  x 

[A2(ti)  sin  /(to)  +  A3(ti)  cos  /(to)]  (4.40) 


4-6 


It  is  now  implied  that  the  const2uits  of  integration  A,,  i  —  1, 2, 3  aire  functions  of  the 
slow  time  ti.  Further,  the  '  now  represents  the  derivative  with  respect  to  the  slow 
time  (ie  ()'  =  dQjdtx).  Integrating  this  equation  yields 


{— A*  +  a  Z2o23o}dto  =  —A\tQ  |q  +  ^ -  6i 


+  C 


\J^(k  -1-  gto) 

1  ^  \ 

+  -^11(^1) 


(4.41) 


62 


T  /  -2  .2  n.  A  A  • 

a  —  lAjCos - Aocos - |-2A2A3sm — | 

9  9  9} 


85 


/IT/  .2  •  Ai  ‘  nA  A 

a,  /  —  1  —  Aj  sin - h  A3  sin  —  +  2A2A3  cos  — 

V  8y  V  9  9  9 ) 


The  Fresnel  Sine  (S)  and  Fresnel  Cosine  (C)  functions  are  as  defined  in  Chapter 
3.  As  can  be  seen,  A'l  is  the  only  secular  term  in  the  solution.  In  order  for  xn  to 
be  bounded  as  t  — »  00,  we  must  have  A\{ti)  =  0.  Further,  application  of  initial 
conditions  for  the  0(1)  system  results  in 


Ai(fi)  =  constant  =  ii(0)  (4.42) 

We  now  “redefine”  Eqs.  4.23  and  4.32. 

xio(to,<i)  =  xi(0)  (4.43) 

k  =  fio-  *3X10  (4.44) 
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The  dependence  of  the  new  “constant”  of  integration,  /lii(ti),  c2ninot  be  determined 
unless  we  go  to  the  0{e^)  system  of  equations.  Finally, 

- — -  Oi  +  C  - — -  62 

\/9  \  I  ^  . 

-  C  b2-S  +  (4.45) 

V9  \  L  . 

The  dependence  of  xn  on  the  fast  time  is  shown  in  the  Fresnel  Sine/ Cosine  terms 
while  its  slow  time  dependence  appears  in  the  constants  of  integration  contained  in 
bi  and  b^. 

We  now  have  a  set  of  coupled  PDF’s  in  the  0(c)  system  of  equations  which 
must  be  solved.  The  equations  are  repeated  here 


®u(fo>^i)  =  S 


dto 

dto 


(*3®10  —  +  *33:11X30  — 


9X20 

dti 


—(*3X10  —  ^)x2i  —  *3X11X20  —  a  X10X20  — 


^X30 

^1 


where 


=  A^{ii)  cos  /(to)  -  A3(<i)  sin  /(<o) 

=  A'2(ti)sin/(to)  +  i4^(ti)cos/(to) 

Cn\ 

This  system  of  equations  can  be  written  as 

(4.46) 

(4.47) 
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Notice  that  at  this  point,  the  B  matrix  is  only  a  function  of  the  fast  time. 


B{to)  = 


0  [13X1(0)  -  fi{to)] 

-  (*3X1(0)  -  fiito)]  0 


(4.48) 


and 


u(fo,<i)  = 


*3X11X30  -  A2  cos  if  {to))  +  A3  sin  (/(to)) 
-*3X11X20  -  oc  X10X20  -  A2sin(/(to))  -  A3  cos  (/(to)) 


(4.49) 


This  system  has  a  solution  of  the  form 


X(to,ti)  =  $(t,0)X(0,0)  +  $(t,0)  jf‘$-'(T)u(T)dr  + 


A2l(ti) 

A3l(ti) 


(4.50) 


The  initial  conditions  for  this  system  will  be  implemented  after  we  complete  the 
derivation  of  the  0(e)  terms.  As  a  result,  we  are  only  interested  in  the  particular 
solution.  A21  and  A31  are  the  constants  of  integration  resulting  from  the  integra¬ 
tion  of  the  system  of  PDE’s.  Because  the  B  matrix  is  only  dependent  on  to,  the 
state  transition  matrix  is  exactly  the  same  as  that  in  the  0(1)  system  of  equations 
(Eq,  4.34).  Multiplication  of  the  inverse  of  the  state  transition  matrix  and  the 
forcing  function  matrix  u  yields  the  following 


<b~^u{to,ti)  = 


-Aj  +  A3i3Xii  -  I  sin  (2/(<o))xioa  A2  +  ^Aoa  x^- 

\  cos  (2/(to))A3a  xio 

-A3  -  A2t3Xii  -I-  5  sin(2/(to))xioa  A3  -  iA2a  xio- 

i  cos  {2f{to))A2Q  Xio 


(4.51) 

In  order  to  achieve  a  uniform  exp2uision,  we  require  that  the  coefficients  of  the  secular 
terms  must  go  to  zero.  The  first,  second  and  fourth  terms  in  each  row  of  the  matrix, 
when  integrated,  will  produce  secular  terms.  The  seculaur  terms  form  the  following 
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system  of  PDE’s 


-i4j  =  -(ia  xio  +  t3^n)A3  (4.52) 

~^3  =  (jQ!  ario  +  *33511)^2  (4.53) 

This  system  allows  us  to  determine  the  dependence  of  Aj  and  A3  on  the  slow  time, 
ti,  however,  the  system  is  extremely  complex.  First,  xn  (Eq.  4.45)  contains  bi  and 
62  which  are  functions  of  Aj  and  A3.  This  makes  the  system  of  equations  nonlinear. 
Second,  xn  also  contains  the  Aii(<i)  term  which  is  unknown  at  this  time  and  can 
only  be  determined  by  solving  the  O(e^)  system  of  equations.  While  the  prospect 
of  deriving  an  approximate  solution  at  this  point  appears  grim,  reasonable  results 
can  still  be  obtadned.  Note  that  the  system  of  equations  for  A2  and  A3  could  easily 
be  solved  if  the  Xn  term  were  a  constant.  In  order  to  make  this  approximation,  we 
must  ignore  all  time  dependence  in  Xn  so  that  it  can  be  decoupled  from  the  xji 
and  X31  system  of  equations.  This  approximation  requires  that  the  variables  A2  and 
A3  given  in  61  and  62  be  given  the  values  of  the  initial  conditions  X2(0)  and  X3(0) 
respectively.  Further,  since  the  time  dependence  in  xn  is  ignored,  the  An  term  is 
also  dropped.  Figure  4.1  depicts  Xn  as  a  function  of  time  for  cin  oblate  gyrostat  with 
a  1  percent  difference  in  transverse  inertias  and  an  initial  state  [0.7,0.4,0.59161,0]. 
The  two  curves  represent  the  gyrostat  with  two  different  non-dimensional  torque 
values  of  0.01  and  0.1.  As  can  be  seen,  the  xn  function  oscillates  and  asymptotically 
approaches  a  constcint  vadue  as  t  — +  00.  By  making  use  of  the  fact  that  both  the 
Fresnel  Sine  and  Fresnel  Cosine  approaich  ^  as  t  —*  00,  this  vedue  can  be  determined 
as 

^llapprox  ~  2^^^  ^ 

where  now 

Cl  =  6i(0)  =  a./^  fx2(0)^cos  ^ — X3(0)*  cos  —  +  2x2(0)x3(0)  sin — | 

V  Sfl'  \  9  9  9 ) 


Ci-C\k 


yi/ir 
'  V9 


(4.54) 
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SW*  Vactor  (.7:.4;.5«161 :0] 


Figure  4.1.  xn  Variation  for  an  oblate  gyrostat  with  initial  conditions  xi  =  0.7, xj  = 
0.4,  X3  =  0.59161  and  e  =  0.01;  g  =  0.01. 


cj  =  6j(0)  =  ^-X2(0)*  sin  ^  +  0:3(0)*  sin  ^  +  2x2(0)x3(0)  cos 


Using  this  approximation  for  Xu  produces  some  excellent  results.  Returning  to  the 
coupled  set  of  PDF’s  given  by  Eqs.  4.52  euid  4.53,  we  have  now  reduced  them  to 
a  constant  coefficient  coupled  set  of  PDF’s  based  on  the  approximation  of  no  slow 
time  dependence  in  Xi.  This  system  of  equations  now  yields  the  following  result. 


■^2(^1)  =  M  cos  (Xti)  +  Az  sin  {Lt\)  (4.55) 

>^3(^1)  =  —A2  sin  {Lti)  +  A3  cos  (Lti)  (4.56) 

where 

L  =  (^a  xio  +  hxuapprox)  (4.57) 
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These  solutions  for  and  A3  as  functions  of  ti  must  now  be  applied  to  the  0(1) 
system  of  equations.  Therefore  Eqs.  4.35  and  4.36  now  become 


a?2o(<o)  =  [^2  cos  {Le  <o)  +  >^3  sin  {Le  to)]  cos  /(<o) 

—  [— sin  {Le  to)  +  A3  cos  {Le  to)]  sin  /(to) 
X3o(to)  =  [A2Cos(£,e  to)  +  A3sin(Ie  to)]sin/(to) 

+  [- A2  sin  {Le  to)  +  A3  cos  {Le  to)]  cos  /(to) 


Further,  with  these  solutions  in  hand,  the  secular  terms  of  the  $“^u(to,ti)  matrix 
(Eq.  4.51)  may  now  be  ignored.  The  integration  from  0  to  t  followed  by  the  multi¬ 
plication  of  the  state  transition  matrix  results  in 


_  (,.53) 
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and 


-  A3C«g  +  tt  +  ^)s[i^] 

+  (7  +  **  +  t) 

_  A2S  j  Sin  (7  I 

The  total  solution  can  now  be  formed. 

Xi  =  xio  +  e  xn  +  0{e^)  (4.60) 

In  expanded  form  this  is 

_ .  f., rv^(*'+9‘)i  . ^\\/^i'‘+st) 

X.  =  x.(0)  +  .|5j^ - - - Jc  +  C^ - - - Jc3 
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The  X2  approximate  solution  is  as  follows. 


Xa  =  ^  ^21  + 


(4.62) 


In  a  slightly  expanded  form  this  is 

xa  =  |i4a  cos  Xi  +  izXnapprox)^  +  ^3  sin  (^-a  xi  +  t3a:ua«>ro*)e  *)]  x 

cos  (^  +  w) 

+  ^i42  sin  xi  +  hxnapprox)^  -  ^Issin  ^(-a  Xi  +  taiiiappro*)^  x 

sin  +  c  xai  +  O(c')  (^-63) 

The  X3  approximate  solution  is  as  follows. 


X3  =  X30  +  £  X31  +  O(e^) 


(4.64) 


=  ^Aa  cos  ^(^a  xi  +  hxnappTox)^  +  -^3  sin  ^2"  "*■  ^)]  ^ 

""  *‘) 


+  j^-Aa  sin  ^(^a  xi  r  isXnapprox)^  +  A3  sin  ^(-a  xi  +  i3Xnappror)£ 


cos 


(t"") 


+  e  X31  +  0{e^) 


(4.65) 


The  solutions  used  for  xai  and  X31  are  not  exact  in  that  we  have  not  solved  for 
Aai(ti)  and  A3i(ti).  This  would  involve  solving  the  C?(e^)  system  of  equations  which 
requires  solving  the  following  integrals  analytically  or  approximately. 


fs{<T) 

Jo 

I  C{<t) 

Jo 


cos  (rdc 
sin  crda 
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Being  unable  to  accomplish  this,  we  assume  this  error  to  be  small  and  proceed. 

To  determine  Aj  and  A3,  we  apply  the  initial  conditions  to  xj  and  X3. 

X2(0,0)  =  X2(0)  =  A2  (4.66) 

X3(0,0)  =  X3(0)  =  A3  (4.67) 

Earlier  we  noted  that  for  the  0(1)  system,  the  norm  of  the  components  of  angular 
momentum  must  equal  one  to  satisfy  the  conservation  of  angular  momentum.  The 
equations  developed  for  iio,  X20,  and  X30  do  satisfy  this  condition.  Since  we  have 
found  Xio  to  be  a  constant,  however,  we  note  that  the  sum  of  squares  of  the  transverse 
angular  momentum  components  must  also  be  constant  for  t  >  0.  This  is  essentially 
the  same  result  that  was  found  by  Sen  and  Badnum  (13).  There  have  been  no 
limitations  on  geometry  to  this  point  and  so  this  result  is  valid  for  near  axisymmetric 
oblate  or  prolate  configurations. 

4.1  Method  of  Multiple  Scales  Results 

We  will  investigate  the  validity  of  this  approximate  solution  for  the  three  types 
of  spinup  problems  described  in  Chapter  2.  Recalling  the  discussion  in  section  2.4, 
we  note  that  there  should  not  be  any  complications  associated  with  instantaneous 
separatrix  crossings  for  an  oblate  gyrostat.  As  seen  in  Figure  2.3,  the  instantaneous 
separatrices  associated  with  unstable  equilibria  migrate  to  the  south  pole  of  the 
momentum  sphere  as  p  increases.  An  oblate  trajectory,  however,  would  begin  and 
stay  in  the  northern  hemisphere  throughout  spinup. 

Dual-spin  gyrostat  dynamics  become  interesting  when  the  initial  condition 
starts  in  the  southern  hemisphere  of  the  momentum  sphere.  A  perfectly  modelled 
system  would  be  able  to  accommodate  the  non-linearity  associated  with  migration  of 
the  instantaneous  separatricqs  in  the  southern  hemisphere.  Recall,  however,  that  by 
choosing  the  difference  in  transverse  inertias  to  be  the  small  parameter,  the  system 
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of  equations  that  describe  the  spinup  are  transformed  from  a  non-linear  system  to  a 
sequence  of  linear  systems.  As  a  restilt,  the  approximate  solution  does  not  contain 
any  of  the  separatrices  associated  with  the  non-linear  system.  Since  both  prolate 
spinup  and  transverse  spinup  have  trajectories  that  interact  with  and  are  influenced 
by  instantaneous  separatrices,  we  expect  to  see  some  degradation  in  the  approximate 
solution  for  these  spinup  problems. 

4.1.1  Oblate  Gyrostats.  Figures  4.2  through  4.7  compare  the  approximate 
solution  versus  the  solution  obt^uned  by  a  4‘^  order  Runge-Kutta  numerical  integra¬ 
tion  scheme.  The  gyrostat  modelled  is  oblate  with  a  1  percent  difference  in  transverse 
inertias  and  an  initial  state  [0.7;  0.4;  0.59161;  0].  The  non-dimensional  torque  value 
is  g  =  0.01.  Again,  fi  is  the  time  scale  for  the  beginning  of  spinup  to  the  end  of 
spinup.  Because  the  scale  of  the  plots  is  large,  the  relative  error  between  each  graph 
is  also  presented.  The  relative  error  is  given  by  the  following  equation. 

Percent  lUbtiv.  Error  =  ~  ^  ,  2  3  (4  jg) 

x,(0) 


The  assumption  of  no  slow  time  dependence  in  the  xi  component  of  angular 
momentum  allowed  the  derivation  of  an  analytic  solution  for  the  terms  of  the  approx¬ 
imate  solution.  Thus,  the  result  for  the  Xi  component  is  identical  to  that  developed 
in  the  straightforward  expansion  method.  Recalling  the  results  from  Chapter  3,  even 
though  this  solution  does  go  in  and  out  of  phase  with  the  numerical  solution,  the 
error  is  less  than  one  percent  due  to  the  small  magnitude  of  the  phase  error  compared 
to  the  initial  value  in  the  xi  component.  Further,  we  can  now  see  that  the  phase 
error  is  caused  by  ignoring  the  slow  time  dependence  in  this  component  of  angular 
momentum.  Figures  4.4  through  4.7  depict  the  X2  and  X3  components  of  angular 
momentum.  Qualitative  results  are  excellent.  Although  there  are  two  curves  in  Fig¬ 
ures  4.4  and  4.6,  only  one  is  visible.  Notice,  however,  the  large  error  at  ~  0.05  in 
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Figiire  4.2.  Numerical  vs  approximate  solution  of  the  Xi  component  of  angular  mo¬ 
mentum  of  an  oblate  gyrostat  with  initial  conditions  xi  =  0.7,  xj  =  0.4, 
X3  =  0.59161  and  e  —  g  —  0.01. 
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Figure  4.4.  Numerical  vs  approximate  solution  of  the  component  of  angular  mo¬ 
mentum  of  an  oblate  gyrostat  with  initial  conditions  xi  =  0.7,  X2  =  0.4, 
X3  =  0.59161  and  e  =  g  =  0.01. 


Figure  4.5.  Percent  relative  error  of  the  X2  component  of  angular  momentum  of  an 
oblate  gyrostat  with  initial  conditions  xi  =  0.7,  X2  =  0.4,  xs  =  0.59161 
and  £  =  g  =  0.01. 
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Figure  4.6.  Numerical  vs  approximate  solution  of  the  X3  component  of  angular  mo¬ 
mentum  of  an  oblate  gyrostat  with  initial  conditions  xi  =  0.7,  X3  =  0.4, 
X3  =  0.59161  and  e  =  g  =  0.01. 
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Figure  4.7.  Percent  relative  error  of  the  X3  component  of  angular  momentum  of  an 
oblate  gyrostat  with  initial  conditions  xi  =  0.7,  X2  =  0.4,  X3  =  0.59161 
and  e  =  g  —  0.01. 
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Figures  4.5  and  4.7  during  the  early  portion  of  the  spinup.  This  “spike”  is  caused  by 
the  assumption  of  a  constant  value  for  the  xn  function  in  the  approximate  solution. 
Further,  the  error  induced  by  the  assumption  manifests  itself  as  phase  error.  To  show 
this,  observe  Figure  4.8  which  depicts  the  amplitude  error  associated  with  the  spinup 
of  this  gyrostat.  We  see  that  the  amplitude  error  is  very  small  in  the  early  portion 
of  the  spinup  maneuver.  Recalling  the  derivation  of  the  approximate  solution,  the 
approximation  made  for  Xn  impacted  the  solution  for  the  slow  time  variation  of  the 
transverse  components  of  angular  momentum.  The  solution  to  Eqs.  4.52  and  4.53  are 
themselves  sines  and  cosines  and  are  the  means  by  which  the  approximate  solution 
captures  the  slowly  varying  frequency  in  the  transverse  components  of  the  system. 
The  amount  of  initial  phase  error  is  related  to  how  far  the  gyrostat  starts  from  the 
oblate  equilibrium  point.  Figure  4.9  depicts  the  Xn  function  for  two  different  oblate 
trajectories.  One  trajectory  begins  very  close  to  the  oblate  equilibrium  point  while 
the  other  begins  very  far  away.  As  seen,  the  approximation  of  the  xu  function  very 
early  in  the  spinup  is  not  as  accurate  as  at  the  end  of  the  spinup.  The  approximation 
gets  worse  as  the  initial  conditions  get  further  away  from  the  equilibrium  point.  The 
same  phenomena  occurs  for  different  values  of  g  (see  Figure  4.1).  A  smaller  non- 
dimensional  torque  value  results  in  a  larger  initial  approximation  error.  Since  this 
error  is  incorporated  as  phase  error  in  the  xs  and  X3  components  of  the  approximate 
solution,  we  can  expect  higher  initial  phase  error  in  the  approximate  solution  for 
oblate  gyrostats  that  begin  spinup  further  away  from  the  oblate  equilibrium  point 
and/or  are  torqued  by  smaller  values  of  g. 

These  results  are  for  one  particular  gyrostat  which  begins  spinup  at  one  initial 
condition.  In  order  to  determine  the  envelope  of  the  solution  for  an  oblate  gyrostat, 
we  look  at  three  different  variables  in  the  equations;  initial  conditions,  torque  (g) 
values  and  differences  in  transverse  inertias  (c). 

4. 1.1.1  Initial  Conditions.  Figure  4.10  depicts  the  maximum  error 
detected  throughout  spin  up  for  a  gyrostat  with  e  =  g  =  0.01  starting  at  various 
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Figure  4.8.  Amplitude  error  for  the  xj  component  of  angular  momentum  of  an 
oblate  gyrostat  with  initial  conditions  xj  =  0.7,  xj  =  0.4,  X3  =  0.59161 
and  £  =  g  =  0.01. 
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Figure  4.10.  Maximum  error  detected  as  a  function  of  initial  conditions  for  an 
oblate  gyrostat. 

initial  conditions  (note  that  the  maximum  error  may  have  occurred  at  different 
times  for  each  component).  The  y  axis  is  the  Xi  initial  condition  with  xj  =  X3  = 
^(1  —  xJ)/2  being  the  transverse  initial  conditions.  Figure  4.10  shows  less  than 
1  percent  error  in  all  three  axes  for  gyrostats  that  start  spinup  from  a  point  half 
way  down  the  north  (oblate)  half  of  the  momentum  sphere.  As  the  initial  condition 
progresses  further  away  from  the  oblate  equilibrium  point,  however,  the  approximate 
solution  rapidly  begins  to  diverge.  This  is  caused  by  the  large  phase  error  induced 
early  in  the  spinup  maneuver  by  assuming  a  constant  value  for  the  xn  function. 
Since  the  primary  interest  in  the  spinup  maneuver  is  the  angular  momentum  at  the 
end  of  spinup,  Figure  4.11  depicts  the  maximum  error  detected  during  the  last  10 
percent  of  the  spinup  maneuver.  As  seen,  with  the  initial  phase  error  absent,  the 
approximate  solution  is  valid  to  within  one  percent  for  the  vast  majority  of  initial 
conditions. 
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Figure  4.11.  Maximum  error  detected  in  the  last  10  percent  of  spinup  as  a  function 
of  initial  conditions  for  an  oblate  gyrostat. 

4. 1.1.8  Torque  Values.  The  torque  value  i^  directly  related  to  the 
validity  of  the  solution  through  Eq.  2.25.  Because  the  approximate  solution  has 
only  incorporated  two  time  scales,  the  solution  is  only  valid  to  t  =  C?(e“^).  The 
reason  is  that  for  t  much  greater  then  C7(e“*),  the  <1  vauriable  ceases  to  be  an  0(1) 
quantity,  and  as  a  result,  the  solution  breaks  down  (11:228*230).  This  fact  limits  the 
torque  values  to  those  which  <ure  larger  than  the  asymmetry  in  the  gyrostat.  The 
unfortunate  result  is  that  this  solution  is  unable  to  model  accurately  the  spinup  of 
gyrostats  with  relatively  large  asymmetry  and  small  spinup  torque.  This  is  not  to 
say,  however,  that  we  cannot  get  good  qualitative  results.  Figures  4.12  and  4.13 
illustrate  this  point.  Here  we  model  the  same  gyrostat  as  in  Figures  4.2  through 
4.7,  but  reduce  the  torque  value  to  jr  =  0.001.  As  can  be  seen,  there  is  excellent 
qualitative  accuracy,  but  a  closeup  reveals  a  slight  phase  error  this  far  into  spinup 
(see  Figure  4.14).  Figure  4.15  illustrates  the  idea  graphically. 

A  second  quzJity  associated  with  the  range  of  vailidity  is  the  error  of  the  ap¬ 
proximate  solution.  In  general,  multiple  scales  solutions  carried  out  to  0{e)  have  an 
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Figure  4.12.  Numerical  vs  approximate  solution  of  the  xj  component  of  angular 
momentum  of  an  oblate  gyrostat  with  initial  conditions  xi  =  0.7, 
X2  =  0.4, 13  =  0.59161  and  e  =  0.01;  5  =  0.001. 

error  of  0(1)  when  t  =  0(€~^).  This  is  because  the  constants  of  integration  which 
develop  in  the  0(1)  terms  are  solved  for  to  0(e)  (cf  Eqs.  4.55  and  4.56).  These 
constants  of  integration  are  of  0(1)  when  t  =  0(e~^)  which  fixes  the  overall  error 
of  the  solution.  Such  an  orderly  generalization  cannot  be  applied  to  the  approxi¬ 
mate  solution  developed  in  this  thesis.  The  approximation  made  for  the  Zn  function 
will  cause  varying  degrees  of  error  as  spinup  time  increases.  Further,  we  employ  an 
approximation  for  the  Fresnel  Sine/Cosine  functions  which  have  their  own  order  of 
error. 


4. 1.1. 3  Difference  in  Transverse  Inertias.  Perturbation  solutions 
assume  that  the  small  parameter,  e,  is  “small.”  While  an  exact  definition  of  “small” 
cannot  be  made,  the  validity  of  the  solution  can  be  demonstrated  for  different  “small” 
values.  Figure  4.16  plots  the  maximum  error  detected  throughout  spinup  for  an 
oblate  gyrostat  with  different  values  of  e.  This  figure  shows  that  we  can  expect  results 
with  less  than  2.4  percent  error  for  an  oblate  gyrostat  that  has  a  5  percent  difference 
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Figure  4.13.  Numerical  vs  approximate  solution  of  the  X3  component  of  angular 
momentum  of  an  oblate  gyrostat  with  initial  conditions  x\  =  0.7, 
X2  =  0.4, 13  =  0.59161  and  e  =  0.01;  jr  =  0.001. 
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Figure  4.14.  Phase  error  that  develops  for  an  oblate  gyrostat  with  g  «  e. 
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Figure  4.15.  Validity  of  solutions  with  g  «  e. 

in  transverse  inertias  and  approximately  1  percent  error  for  a  gyrostat  that  has  a 
2  percent  difference  in  transverse  inertias.  The  data  in  this  figure  must  be  viewed 
carefully,  however,  since  this  gyrostat  has  started  with  an  initial  condition  fairly  close 
to  the  oblate  equilibrium  point.  We  know  the  error  values  will  increase  as  the  initial 
conditions  move  further  from  the  oblate  equilibrium  point.  In  contrast,  recall  that 
the  maximum  error  associated  with  oblate  gyrostats  stating  spinup  further  away 
from  the  oblate  equilibrium  point  is  caused  by  the  phase  error  induced  early  in  the 
spinup  maneuver.  The  error  is  considerably  less  when  the  comparisons  are  confined 
to  the  last  10  percent  of  the  spinup  maneuver. 

4. 1.1.4  Cone  Angle  Approximations.  As  discussed  in  Chapter  2, 
the  cone  (or  nutation)  angle  is  the  angle  formed  by  the  total  angular  momentum 
vector  and  the  symmetry  axis  of  the  gyrostat.  The  equation  for  this  angle  is  given 
by  Eq.  2.23.  Throughout  the  analysis  for  oblate  gyrostats,  the  maximum  error 
associated  with  the  Xi  component  of  angular  momentum  has  been  less  then  one 
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Figure  4.16.  Maximum  error  detected  as  a  function  of  s  for  an  oblate  gyrostat  with 
initial  conditions  xi  =  0J,X2  =  0.4,  X3  =  0.59161  ;</  =  e. 

percent.  From  Figure  4.10,  this  maximum  error  occurred  for  an  oblate  gyrostat 
with  a  one  percent  difference  in  transverse  inertias  but  starting  just  north  of  the 
equator  of  the  momentum  sphere.  From  Figure  4.16,  this  maximum  error  occurred 
for  an  oblate  gyrostat  starting  fairly  close  to  the  oblate  equilibrium  point  but  with 
a  ten  percent  difference  in  transverse  inertias.  The  essential  result,  however,  is 
that  Eq.  4.61  is  a  reasonable  approximation  for  near  axisymmetric  oblate  gyrostats 
which  begin  the  spinup  maneuver  with  initial  conditions  in  the  northern  half  of  the 
momentum  sphere.  Further,  this  equation  can  be  used  to  determine  the  cone  angle 
at  any  time  during  the  spinup  maneuver. 

4.1.2  Prolate  Gyrostats.  As  mentioned  earlier,  our  selection  of  the  differ¬ 
ence  in  transverse  inertias  allows  us  to  solve  linear  systems  of  equations  for  the  terms 
of  the  approximate  solution.  Because  the  systems  of  equations  are  linear,  our  ap¬ 
proximate  solution  does  not  incorporate  any  of  the  separatrices  associated  with  the 
original  nonlinear  governing  equations.  The  result  is  that  a  phase  error  will  develop 
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as  the  instantaneous  separatrices  begin  to  migrate  tow2Lr<l  the  prolate  equilibrium 
point. 

To  see  this,  recall  that  in  a  ^  =  0  system,  the  separatrix  separates  different 
kinds  of  motion.  In  Figure  4.17,  we  present  a  cartoon  of  this  concept  in  the  xj,  X3 
pl2me.  It  is  clear  that  the  phase  of  the  trajectory  inside  the  separatrix  is  different  than 
that  on  the  outside.  As  /i  increases,  the  prolate  equilibrium  point  will  go  through  two 
bifurcations  as  discussed  in  Chapter  2,  and  the  phase  of  the  trajectories  on  either  side 
of  the  separatrix  will  change.  The  net  result  is  that  the  instantaneous  separatrices 
in  the  nonlinear  system  influence  the  phase  of  the  trajectories  on  either  side  of  the 
separatrix.  Since  the  approximate  solution  does  not  contain  any  sep2u^atrices  to 
influence  the  trajectories,  we  expect  a  phase  error  to  develop.  We  also  expect  the 
magnitude  of  the  phase  error  to  be  larger  when  the  separatrix  is  near  the  equator  of 
the  momentum  sphere  since  the  trajectory  paths  are  longer. 

Unfortunately,  the  phase  error  induced  by  the  instantaneous  separatrices  is 
not  the  only  phase  error  in  the  approximate  solution.  Just  as  for  an  oblate  gyrostat, 
the  Xii  approximation  induces  an  additional  phase  error.  The  combination  of  these 
errors  greatly  hinders  the  usefulness  of  this  approximate  solution  for  prolate  spinup. 

Figures  4.18  through  4.20  illustrate  the  spinup  of  a  prolate  gyrostat.  We  can 
see  that  the  approximate  solution  begins  to  get  out  of  phase  near  /i  =  0.1.  Un¬ 
fortunately,  because  the  approximate  solution  begins  with  a  relatively  large  error, 
which  is  compounded  by  an  additional  phase  error  from  its  inability  to  model  the 
instantaneous  separatrices,  the  remainder  of  the  solution  is  flawed.  Notice  that  the 
gyrostat  begins  spinup  faurly  close  to  the  prolate  equilibrium  point  {P^).  Because 
the  Xii  approximation  gets  worse  as  the  initial  conditions  move  further  from  P^,  we 
do  not  expect  any  results  better  than  those  shown. 

As  discussed  in  Section  2.4,  the  last  bifurcation  associated  with  the  instanta¬ 
neous  separatrices  occurs  at  ft  equal  to  the  larger  of  |i2|  or  liaj.  After  this  point,  the 
momentum  sphere  contains  only  two  equilibrium  points,  both  of  which  axe  centers. 
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Figtire  4.17.  Cartoon  of  different  trajectory  paths  in  a  ^  =  0  nonline2ir  system. 


Figure  4.18.  Numerical  vs  approximate  solution  of  the  Xi  component  of  angular 
momentum  of  a  prolate  gyrostat  with  initial  conditions  xi  =  0.9,  xs  = 
0.30822,  X3  =  0.30822  and  c  =  ^  =  0.01. 
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Figure  4.19.  Numerical  vs  approximate  solution  of  the  Xj  component  of  angular 
momentum  of  a  prolate  gyrostat  with  initial  conditions  xi  =  0.9,  X2  = 
0.30822,  X3  =  0.30822  and  e  =  i7  =  0.01. 
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Because  there  are  no  more  instantaneous  separatrices  to  influence  the  trajectories, 
the  approximate  solution  should  well  approximate  the  system  again  after  the  last 
bifurcation  occurs.  Figures  4.21  through  4.23  illustrate  this  idea.  These  graphs 
represent  the  error  between  the  two  solutions.  We  used  the  numerical  results  at 
H  =  0.505  as  the  initial  conditions  for  the  last  portion  of  the  spinup.  These  graphs 
also  clearly  show  the  early  phase  error  associated  with  the  xn  approximation  and 
the  instantaneous  separatrices.  The  small  error  seen  during  the  last  portion  of  the 
spinup  is  on  the  same  order  as  that  seen  during  oblate  spinup.  Recognize  that  if  we 
use  the  numerical  results  anywhere  from  when  the  solution  starts  to  get  out  of  phase 
to  fi  —  the  larger  of  (tsl  or  {tal  as  initial  conditions,  the  approximate  solution  will  not 
be  as  accurate.  The  reason  is  that  even  though  the  approximate  solution  is  starting 
from  a  correct  initial  condition,  it  does  not  have  the  correct  frequency  information 
until  the  effects  of  the  instantaneous  separatrices  are  eliminated. 

4.1.3  Transverse  Spinup.  Transverse  spinup  should  experience  the  same 
problems  associated  with  prolate  spinup  except  the  errors  will  be  more  pronounced. 
Recall  that  in  this  case,  the  trajectory  starts  near  a  stable  flat  spin  equilibrium  point. 
Further,  unless  the  torque  vaJue  is  extremely  small,  the  trajectory  will  have  too  much 
energy  and  cross  an  instantaneous  separatrix  almost  immediately.  This  will  result 
in  a  large  phase  error  between  the  approximate  and  numerical  solutions.  The  error 
is  the  composite  of  the  Xn  approximation,  which  is  the  worst  at  small  /z,  and  the 
magnitude  of  the  phase  error  which  is  largest  near  the  equator.  Figures  4.24  through 
4.26  compare  the  approximate  2uid  numerical  solutions  for  a  prolate  gyrostat  starting 
at  the  exact  transverse  equilibrium  point.  Figures  4.27  through  4.29  show  the  same 
gyrostat  starting  spinup  just  offset  from  the  equilibrium  point.  In  each  case,  the 
instantaneous  separatrix  crossing  occurs  almost  immediately.  This  is  evidenced  by 
the  fact  that  the  Xj  and  X3  components  immediately  begin  to  oscillate  between  -1 
and  1. 
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Figiire  4.20.  Numerical  vs  approximate  solution  of  the  X3  component  of  angular 
momentum  of  a  prolate  gyrostat  with  initial  conditions  xi  =  0.9,  X3  = 
0.30822,  X3  =  0.30822  and  c  =  9  =  0.01. 


sum  \MClor-(.0:.30822;.30e22;0] 


Figure  4.21. 


Percent  relative  error  in  the  xi  component  of  angular  momentum  of 
a  prolate  gyrostat.  Initial  conditions  are  matched  at  /x  =  0  2uid  fi  = 


0.505. 
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Figure  4.22. 


Percent  relative  error  in  the  Xj  component  of  angular  momentum  of 
a  prolate  gyrostat.  Initial  conditions  are  matched  at  p  =  0  and  p  = 


0.505. 


Stua  V«CtOr«4J»;.30822-..30e22»] 


Figure  4.23. 


Percent  relative  error  in  the  X3  component  of  angular  momentum  of 
a  prolate  gyrostat.  Initial  conditions  are  matched  at  /x  =  0  and  p  = 


0.505. 
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Figure  4.24.  Numeric2kl  vs  approximate  solution  of  the  xi  component  of  angular 
momentum  of  a  prolate  gyrostat  during  transverse  spinup. 

Recall  that  our  solution  is  only  valid  for  g  >  e.  In  order  to  make  use  of 
this  approximate  solution,  we  would  need  to  model  extremely  near  axisymmetric 
gyrostats  {e  =  0.00001)  with  extremely  small  torque  values.  This  is  not  very  practical 
since,  in  the  limit,  it  approaches  the  exact  solution  for  e  =  ^  =  0  (i.e.  axisynunetric 
and  no  spinup). 

4.2  Summary 

The  approximate  solution  developed  in  this  chapter  essentially  employs  the 
method  of  multiple  scales.  The  solutions  are  not  complete  in  that  we  are  unable 
to  solve  the  Oie^)  system  of  equations  to  determine  the  constants  of  integration  in 
the  0(e)  system.  Further,  in  order  to  derive  an  analytical  expression  for  the  0(1) 
system,  we  ignore  the  slow  time  dependence  in  the  Xi  component  of  angulau’  momen¬ 
tum.  While  this  approximation  leads  to  large  phase  errors  in  Xi,  the  magnitude  of 
the  error  is  small  compared  to  the  magnitude  of  the  overall  component  of  angular 
momentum.  The  result  is  an  overall  small  error  in  the  Xj  component.  The  benefit 
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Figure  4.25. 


Figure  4.26. 


MMa  v«cior>(0;0;l  :0I 


Numerical  vs  approximate  solution  of  the  component  of  angular 
momentum  of  a  prolate  gyrostat  during  transverse  spinup. 


MM*  v«ctor-(0;0;l  ;0) 


Numerical  vs  approximate  solution  of  the  X3  component  of  angular 
momentum  of  a  prolate  gyrostat  during  tr2Lnsverse  spinup. 
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Figt  re  4.27.  Numerical  vs  approximate  solution  of  the  xi  component  of  angular 
momentum  of  a  prolate  gyrostat  during  transverse  spinup. 


atat*  vw;lorM.30822;.30822;.9;0] 


mu 


Figure  4.28.  Numerical  vs  approximate  solution  of  the  X2  component  of  angular 
momentum  of  a  prolate  gyrostat  during  transverse  spinup. 
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Figure  4.29.  Numerical  vs  approximate  solution  of  the  X3  component  of  angular 
momentum  of  a  prolate  gyrostat  during  transverse  spinup. 

of  the  approximation  is  that  it  allows  the  analytical  development  of  the  transverse 
components  of  angular  momentum. 

We  compared  the  approximate  solution  to  the  three  spinup  problems  discussed 
in  Chapter  2.  For  oblate  spinup,  the  approximate  solution  provides  excellent  results. 
For  oblate  gyrostats  that  begin  spinup  close  to  the  oblate  equilibrium  point,  the  error 
in  all  three  components  of  angular  momentum  are  very  small.  The  reason  is  that  the 
influence  of  the  instantaneous  separatrices  are  initially  small  and  remain  small  as  the 
spinup  maneuver  continues.  As  the  initial  conditions  for  the  gyrostat  get  close  to 
the  equator  of  the  momentum  sphere,  the  influence  of  the  instantaneous  sepauratrices 
are  initially  felt  and  result  in  a  large  initial  phase  error.  However,  as  the  spinup 
maneuver  progresses,  the  instantaneous  separatrices  migrate  further  away  from  the 
trajectory  and  the  error  values  decrease.  Additionadly,  we  discussed  the  validity  of 
the  approximate  solution  in  terms  of  time  amd/or  torque  vadue.  Because  we  only 
use  two  time  scales,  the  approximate  solution  is  only  valid  to  <  =  This 

implies  that  the  torque  must  be  greater  tham  the  e  value  to  achieve  quantitatively 
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accurate  solutions.  However,  qualitatively  accurate  solutions  are  available  out  to  t  = 
Finally,  we  determined  that  oblate  gyrostats  with  a  one  percent  difference 
in  transverse  inertias  are  well  approximated. 

For  prolate  spinup,  the  approximate  solution  is  hindered  by  phase  errors  arising 
from  the  in  approximation  and  its  inability  to  account  for  instantaneous  separatri- 
ces.  This  results  in  an  “instantaneous  separatrix  window”  in  which  the  approximate 
solution  is  not  able  to  compensate  for  the  frequency  change  as  much  as  the  nonlinear 
system.  The  “window”  closes  after  the  last  bifurcation  occurs  at  the  south  pole  of 
the  momentum  sphere.  Numerical  results  for  the  components  of  angular  momentum 
at  any  time  after  the  last  bifurcation  are  valid  initial  conditions  to  approximate  the 
remainder  of  the  spinup. 

The  approximate  solution  is  not  well  suited  to  model  transverse  spinup.  The 
torque  values  required  foi  accurate  use  of  the  approximate  solution  are  so  small,  that 
the  solution  is  rendered  impractical. 
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V.  Cylindrical  Coordinates 

The  approximate  solution  derived  in  Chapter  4  provides  excellent  results  for 
oblate  gyrostats  with  initial  conditions  that  start  in  the  upper  half  of  the  momentum 
sphere.  The  approximate  solution  begins  to  breakdown,  however,  due  to  the  phaae 
error  introduced  by  assuming  a  constant  value  of  the  xn  function.  This  is  most 
clearly  demonstrated  in  Figure  4.10.  The  maximum  error  that  begins  to  develop 
in  the  transverse  components  of  angular  momentum  is  due  to  this  approximation. 
We  have  noted  that  this  error  occurs  early  in  the  spinup  maneuver  where  the  Xn 
approximation  is  least  valid.  The  idea  of  this  chapter  is  to  derive  an  approximate 
solution  that  does  not  suffer  this  deficiency. 

A  different  approximate  solution  was  developed  by  Lt  Col  William  P.  Baker, 
Associate  Professor  of  Mathematics,  Air  Force  Institute  of  Technology  and  is  included 
here  for  evaluation  and  completeness.  Lt  Col  Baker  observed  that  the  transverse  an¬ 
gular  momentum  components  are  characterized  by  constant  amplitude  but  varying 
frequency.  His  idea  was  to  transform  the  system  of  equations  to  cylindrical  coor¬ 
dinates  and  then  use  the  method  of  multiple  scales  to  derive  another  approximate 
solution. 

To  develop  the  solution  we  proceed  as  follows.  Eqs.  2.18  through  2.20  with  the 
substitution  of  as  =  —  *3  reduce  to 


Xi 

=  aex2X3 

(5.1) 

X2 

=  (*3Xl  -  H)X3 

(5.2) 

X3 

=  -(*3Xi  -  n)X2  -  asxiX2 

(5.3) 

A 

=  9 

(5.4) 
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Let 


Xi 

=  Xi 

(5.5) 

X2 

=  pcos  ^ 

(5.6) 

X3 

=  psin^ 

(5.7) 

where  p  is  the  amplitude  and  <l>  is  the  phase  auigle  of  the  system.  The  system  is 
subject  to  the  following  initial  conditions 

xi(0)  =  xi(0) 


<>(0)  =  +  4(0) 


To  get  the  system  of  equations  into  a  form  from  which  we  can  employ  a  multiple 
scales  approach,  we  start  with 


p*  =  4  +  4 

(f>  =  tan"^(x3/x2) 

The  time  derivative  of  these  expressions  are 


X2X2  +  X3X3 

— QexiX2X3  =  — aexip*cos<^sin^ 

1  X3X2  -  X2X3 

1  +  (X3/X2)^  x| 

One  further  manipulation  is 

=  X3X2  —  X2X3  =  — (bxi  —  p)p^  —  aexip^  cos^  <l> 


pp  = 

^  = 


(5.8) 

(5.9) 
(5.10) 


(5.11) 
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Substituting  Eqs.  5.6  and  5.7  into  Eq.  5.1  and  using  Eq.  5.9  and  5.11  results  in  the 
following  system  of  equations. 


Xl 

P 


2  ...  2  .  I 

aep  cos  9  sm  <p=^-prp  sm  2<p 

05 

— Qexipcos^sin^  =  — — xipsin2^ 


X  <3(e  ae 

=  -(*3X1  -  p)  -  yXi  -  yXi  COS 


We  simplify  these  equations  further  by  letting 


^  =  2<l> 

-  Ct£ 

and  obtain  the  following  system  of  equations. 

Xl  =  6p^s\ntl> 
p  =  — ^xipsin^ 

tp  =  — 2(*3Xi  —  p)  —  2^Xi  —  26xi  cos  0 
We  seek  solutions  of  the  form 


(5.12) 

(5.13) 


(5.14) 

(5.15) 

(5.16) 


xi(to,  ti;^)  =  (rio(toi^i)  +  ^(Cii(tojti)  +  ^*2^12(^05^1)  +  (5*17) 

p(to,  ti;^)  =  Po(to5^i)  +  ^Pi(<05ti)  +  ^*P2(^ojti)  +  0(^^)  (5.18) 

t(^{tQ]6)  =  i^o(to)  +  5*^i(to)  +  ^*^2(^0)  +  0{6^)  (5.19) 

Again  the  fast  time  “<0”  is  defined  as  t  while  the  slow  time  “ti”  is  defined  as  6t.  Note 
that  we  msdce  the  initial  assumption  of  no  slow  time  dependence  in  the  phase  angle. 
This  assumption  will  be  discussed  in  the  final  form  of  the  approximate  solution.  We 
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employ  the  following  differential  operator  to  incorporate  the  different  time  scales. 

±_d_  ^ 

dt  dto  ^  dt\ 

Further,  we  still  must  conserve  angul2ir  momentum  and  so  we  impose  the  constraint 

||xll^  =  l,Vt  >0  (5.21) 


This  results  in  the  same  orthogonality  condition  on  adl  0(e)  terms  as  seen  in  the 
cartesian  solution.  Substituting  Eqs.  5.17  through  5.20  into  Eqs.  5.14  through  5.16 
and  collecting  0(1)  terms  yields 

0  (5.22) 

0  (5.23) 

-2(*3Xio  -  m)  (5.24) 

This  system  is  subject  to  the  following  initial  conditions. 


ox  10 
dto 
dpo 
dto 
d^Q 
dto 


a;io(0,0)  =  xi(0) 

po(0,0)  =  p{0) 

MO)  =  2m 
||xop  =  l,v  t  >  0 


Note  that  the  constraint  refers  to  the  total  component  of  angular  momentum.  The 
solution  to  the  0(1)  system  with  the  application  of  initial  conditions  is 


Xio(to,ti)  =  v4o(fi)  =»  xio(0,0)  =  Ao(0)  =  xi(0) 
Po{to,h)  =  />o(^i)  Po(0,0)  = /)o(0)  =  p(0) 
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(5.25) 


rh 

rpoito)  -  /  [2y<7  +  2/io  -  213X10]  da 

Jo 

=  gtl  +  2  (/lo  —  *3X10)  to  +  $0 

where  60  is  defined  so  that  the  initial  conditions  are  satisfied. 

60  =  2tan~^  (x3(0)/x2(0)) 

To  simplify  the  equations,  we  make  the  following  definitions. 

k  =  /^o- *3X1(0) 

/(to)  =  +  2kfo  +  ^0  V’o(to)  = /(to) 


The  0(S)  system  of  equations  is 

_ 

Oto 

dpi  _ 
dto 

d^i  _ 
dto 


2  .  .  ^10 

p^smipo  - 

C/t\ 

.  ,  dpo 

-Xiopo8mV»o  -  -ST 

0\t\ 

-2t3Xii  —  2xio  -  2xio  cos  ^o 


subject  to  the  following  initial  conditions 


xu(0,0)  =  0 

Pi  (0,0)  =  0 

V'i(O)  =  0 

xo  •  xi  =  0,  V  f  >  0 


We  solve  Eq.  5.29  to  find  the  slow  time  dependence  of  Aq. 

xii  =  [po(ti)^  sin  (ga^  +  2ka  +  9o)  -  /lo'(ti)|  da 


(5.26) 


(5.27) 

(5.28) 


(5.29) 

(5.30) 

(5.31) 


(5.32) 
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=  Poihy  1^  jcos  (k^/g  -  0o)s 

-  C  +  5^0)  j  sin  (le^lg  -  Oq) 

-  co8(iV5-^o)5 

+  C  ^-^(A:)^  sin  (^k^/g  —  ^o)H  ~  •^*(^1)^0  +  -^iC^i) 


(5.33) 


The  “constant”  of  integration  for  this  equation  is  /li(ti).  Again,  the  term  “constant” 
indicates  that  the  term  is  constant  with  respect  to  to-  The  complete  form  of  the 
solution  is  as  given,  however,  to  simplify  the  remainder  of  the  analysis  we  use  the 
following  definitions. 


^oif)  =  /  coafda 

Jo 

=  ^  jcos  (^k^/g  -  do)c 
+  5  sin  -  ^0) 

-  cos  {k^jg  - 

-  5  sin  -  ^0)  I  (5.34) 

(5.35) 

So{f)  =  /  sin/dt 

Jo 

=  ^  jcos  {Je^fg  -  0o)  S  +  fl'^o) j 

-  C  +  ^^0)^  sin  (fc^/sr  -  ^0) 

-  cos  (^k'^/g  -  do)S 

+  C  sin  {y/g  —  0o)|  (5.36) 
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where 


/(<t)  =  +  2ka  +  (5.37) 

The  Fresnel  Sine  and  Fresnel  Cosine  functions  are  as  defined  in  Chapter  3,  and 
their  integrals  are  developed  in  detail  in  Appendix  B.  Eaurh  term  in  the  xn  solution 
is  bounded  with  the  exception  of  In  order  to  have  a  completely  bounded 

solution,  must  equal  zero  and  by  application  of  the  initial  condition 

Ao(ti)  =  constant  =  Xi(0)  (5.38) 


The  resulting  solution  for  xn  is 


Xii(to5^i)  =  Po{h)So{f)  +  Ai(ti) 

(5.39) 

We  now  solve  Eq.  5.30. 

Pi  =  jf  [-xi(0)po(<i)  sin  (gfcr^  +  2k<T  +  9o)  -  Po'(<i)] 

(5.40) 

This  equation  is  very  similar  to  Eq.  5.32.  Again,  in  order  to  have  bounded  solutions, 
we  must  have  Po(fi)  =  0  and  by  application  of  the  initial  condition 

Po(^i)  =  constant  =  p(0) 

(5.41) 

As  a  result 

Pi(foj^i)  =  ~xi(0)p^(0)5o(/)  + 

(5.42) 

Turning  our  attention  to  Eq.  5.31,  we  write  it  in  terms  of  the  Eqs.  5.39  and 
5.25  which  have  been  derived. 


dipt 

dto 


-2z3P^(0)5o(/)  -  2t3Ai(fi)  -  2xi(0)  -  2xi[0)cosf 


(5.43) 
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We  employ  the  following  relation  (see  Appendix  B). 


/  So  if  (a))  da 
Jo 

=  tSo  (/(f))  +  ;r  (cos  /  -  cos  9o)  +  -So 


Thus,  the  solution  is 


V'i(fo)  =  -2t3P^(0)  fo5o(/)  +  ^(cos/-cosdo)  + j5o(/) 
—  2fo  (*3Ai(fi)  +  xi(0)]  —  2xi(0)Co(/)  + 


(5.44) 


The  constant  of  integration  is  ^  and  by  application  of  the  initial  condition  we  see 
that  /d  =  0.  With  the  0{6)  system  of  equations  solved,  we  must  determine  the  nature 
of  the  constants  of  integrations  Ai(fi)and  pi(fi).  Note  that  from  the  orthogonality 
condition  specified  on  the  0(6)  system,  we  require  that 


Xo  •  Xi  =  0,  V  f  >  0 


Performing  this  product  results  in 


xi(0)Ai(fi)  +  p(0)pi(fi)  =  0 


(5.45) 


Therefore,  if  Ai(fi)  =  0,  we  must  require  pi(fi)  =  0.  In  order  to  determine  Ai(fi), 
we  proceed  to  the  system  of  equations. 


^12 

^0 

dpi 

dtf) 

drjfi 

dto 


2popi  sin  ipo  +  pQ^\  cos 

-(xioPi  +  xupo)  sin  V’o  -  PoXioV’i  cos - 


dpi 

di\ 


-2*3X12  —  2xio  -  2xio  cos  V’o  -  2x1001  cos  0o  -  2xio  sin  0o 


(5.46) 

(5.47) 

(5.48) 
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This  system  is  subject  to  the  following  initial  conditions. 


xn(0,0) 

/»^(0,0) 

^^>2(0) 

||xi1|^  +  2xo-X2 


=  0 
=  0 
=  0 

=  0,  V  <>  0 


Expressing  Eq.  5.46  in  terms  of  known  functions 


dx 


12 


3to 


-2p^(0)xi(0)5o(/)  sin  /  +  2p(0)pi(<i)  sin  / 
2i3p"(0) 

2p^(0)  [ia^iCti)  +  xi(0)]  tocosf 
2xx(0y(0)Co(/)cos/-M<i) 


toSoif)  +  ^  (cos  /  -  cos  do)  +  -Soil) 
^9  9 


cos  / 


The  solution  to  Eq.  5.46  is 


xi2{<o,ti)  =  -p'(0)xx(0)[5„(/)]'  +  2p(0)pi(tx)5„(/) 


25 

2V(0),. 


(*3>li(/i)  +  Xx(0))  [sin  /  -  sin  do] 
(*3^l(tl)  +  Xi(0))  Co{f) 


2  .  2i3p'(0) 


25 


-  Xx(0)p^(0)[a(/)]^  + 

+  +  -^2(^1) 


cosdoCo(/) 


(5.49) 


(5.50) 


Every  term  in  the  solution  is  bounded  except  i4j(ti).  Therefore  we  require  that  it 
be  equal  to  zero  and  by  application  of  initial  conditions  we  find  that  Ax(<x)  =  0. 


5-9 


Further,  by  Eq.  5.45  we  find  that  pi(fi)  must  also  equal  zero.  However,  we  show  the 
same  process  that  we  applied  to  Eq.  5.47  to  verify  this. 


=  -poiPo  -  xi^)5o(/)  -  Xipi  sin  / 

Cno 

r  1  k 

+  2i3Po^i  toSoif)  +  —  (cos  /  -  cos  9o)  +  -So{f)  cos  / 

I  Ig  g 

+  2xi*/»o^o  cos  /  +  2xi*poCo(/)  cos  /  -  (5.51) 

Using  the  same  relations  that  solved  Eq.  5.46  we  have  the  solution 

(.)  =  i^O)  [/(O)  -  x,'(0)]  |S4/)1^  -  X.(0)^.(«.)S.(/) 

+  2V(0)x,(0)  ^S.(/)sm/+^C<,(2/) 

-  ^  cos  «oC„(/)  +  [sin  f  _  sin  «„) 

-  H-  [Cs(/)1- 

g  i 

-  toPif{ti)  +  p^ih)  (5.52) 

This  solution  is  the  same  as  that  for  the  Xi2  solution  in  that  all  terms  are  bounded 
except  p'x{ti).  Therefore  pi(fi)  =  0  amd  by  application  of  initiad  conditions,  we  find 
that  pi(ti)  =  0  a«  required.  With  the  constamts  of  integration  determined,  we  now 
return  to  Eqs.  5.39  and  5.42  aind  note  their  finad  solution. 

Xu(to,ti)  =-pimSoU)  (5.53) 

fh{U,tx)  = -xi(0)po(0)5„(/)  (5.54) 

To  summatrize,  the  approximate  solution  to  0{e)  using  cylindricad  coordinates 
is  as  follows. 

Xl(t)  =  X,  -  ypg(0)So(/)  +  0(£»)  (5.55) 
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(5.56) 


fit)  =  «)(0)  -  y«.(0)MO)S<,(/)  +  0(£’) 

^(0  =  ^  ^  ~  |*3p^(0)  ^tSoif)  +  ^  (cos  /  -  cos  ^o) 

+  j5o(/)j  +  /x,(0)  -  xi(0)C7<,(/)||  +  0(6^)  (5.57) 

Finally 


X2(t)  =  p(t)  cos  </)(t)  (5.58) 

X3(t)  =  p(t)  sin  </>(t)  (5.59) 


Recognize  that  to  0(1),  we  have  satisfied  the  conservation  of  angular  momen¬ 
tum  requirement.  In  the  beginning  of  the  analysis,  we  made  the  assumption  of  no 
slow  time  dependence  in  the  phase  angle.  Careful  analysis  of  Ek}s.  5.55  through  5.57 
reveal  some  interesting  results.  Although  the  equations  for  Xi  and  p  were  derived 
with  the  assumption  of  some  slow  time  dependence,  it  never  shows  up  in  the  fi¬ 
nal  form  of  the  solution.  However,  slow  time  dependence  is  implicitly  found  in  the 
equation  for  (f>  (Fq.  5.57)  even  though  it  was  not  assumed  to  be  there  initially.  As 
a  result,  by  using  cylindrical  coordinates,  the  approximate  solution  is  essentially  a 
series  expansion  of  the  system  of  equations  using  a  single  time  scale. 

5.1  A  Comparison  Between  the  Cartesian  and  Cylindrical  Approximate  Solutions 

5.1.1  Oblate  Gyrostats.  Figures  5.1  through  5.3  present  the  relative  error 
associated  with  each  component  of  anguW  momentum  for  both  the  cartesism  and 
cylindrical  solutions.  The  model  is  a  near  axisymmetric  (e  =  0.01)  oblate  gyrostat 
with  initial  state  [0.7;  0.4. ;  0.59161;  0].  As  seen,  the  solutions  for  the  xi  component  are 
identical.  This  is  expected  since  the  underlying  assumptions  between  the  solutions 
are  the  same.  Specifically,  the  initial  conditions  for  Xj  and  X3  are  used  in  both 
forms.  Further,  both  solutions  use  the  assumption  of  no  slow  time  dependence.  The 
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transverse  components  show  a  dramatic  difference  in  the  early  portion  of  the  spinup 
maneuver.  The  cylindrical  solution  does  not  exhibit  the  “spike”  that  the  cartesian 
solution  displays.  In  order  to  examine  this  finding  in  more  detail,  we  compare  the 
two  solutions  using  the  same  oblate  gyrostat  starting  at  a  point  further  from  the 
oblate  equilibrium  point  (Figures  5.4  and  5.5).  Here  we  observe  the  same  behavior 
early  in  the  spinup,  but  also  notice  that  the  cartesian  solution  displays  less  error 
than  the  cylindrical  solution  towards  the  end  of  the  spinup  maneuver.  To  explain 
this  behavior,  recall  the  development  of  each  solution.  In  the  cartesian  solution,  the 
“spike”  is  the  result  of  phase  error  caused  by  the  approximation  of  a  constant  value 
for  the  Xii  function.  The  cylindrical  solutions  makes  no  such  approximation  and  so 
there  should  not  be  any  significant  phase  error  aissociated  with  this  region  of  the 
spinup.  Towards  the  end  of  the  spinup  maneuver,  the  lower  error  of  the  C2u^tesisui 
solution  is  the  result  of  the  Xn  approximation.  This  illustrates  an  interesting  result 
between  the  two  solutions.  As  noted  in  Chapter  4,  the  Xn  approximation  impacts 
the  transverse  components  in  that  it  captures  the  slowly  varying  frequency  of  the 
system.  The  cartesian  solution  is  based  on  a  fast  and  a  slow  time  scale.  The  slow 
time  scale  captures  the  phase  or  frequency  of  the  system.  By  assuming  the  const^ult 
value  for  the  Xn  function  was  Xiioppro»(f  =  oo)j  the  frequency  axljustment  of  the 
system  is  incorporated  in  the  transverse  components  for  t  — ♦  oo.  As  seen,  this 
phase  correction  is  not  as  accurate  in  the  early  portion  of  the  spinup,  but  increases 
in  accuracy  as  time  increases.  The  cylindrical  solution  implicitly  incorporates  slow 
time  dependence  but  it  does  not  benefit  from  the  t  —*  oo  phase  correction  that  the 
cartesian  solution  has.  As  a  result,  the  cylindrical  solution  has  a  larger  error  in  the 
latter  portion  of  the  spinup  maneuver. 

To  illustrate  this  point  in  a  different  light,  observe  Figure  5.6.  Here  we  compare 
the  maximum  error  detected  between  cartesiain  and  cylindrical  solutions  in  the  last 
ten  percent  of  spinup  for  a  near  axisymmetric  (e  =  0.01)  oblate  gyrostat  starting  at 
incremental  points  on  the  momentum  sphere.  The  initial  Xi  component  is  the  y  axis 
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Figure  5.2.  Comparison  between  cartesian  and  cylindrical  equations  for  percent 
relative  error  of  the  X2  component  of  angular  momentum  of  an  oblate 
gyrostat  with  initial  conditions  xi  =  0.7,  X2  =  0.4,  X3  =  0.59161  and  e  = 

g  =  0.01. 
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Figure  5.3. 


Figure  5.4. 
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Comparison  between  cartesian  and  cylindrical  equations  for  percent 
relative  error  of  the  X3  component  of  angular  momentum  of  an  oblate 
gyrostat  with  initial  conditions  Xi  =  0.7,  Xj  =  0.4,  X3  =  0.59161  and  c  = 

g  =  0.01. 


Comparison  between  cartesian  and  cylindrical  equations  for  percent  rel¬ 
ative  error  of  the  X2  component  of  angular  momentum  of  an  oblate  gyro¬ 
stat  with  initial  conditions  xi  =  0.4,  X2  =  0.64807,  X3  =  0.64807  and  e  = 

g  =  0.01. 
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Figure  5.5.  Comparison  between  cartesian  and  cylindrical  equations  for  percent  rel¬ 
ative  error  of  the  X3  component  of  angular  momentum  of  an  oblate  gyro¬ 
stat  with  initial  conditions  xi  =  0.4,  Xs  =  0.64807,  X3  =  0.64807  and  e  = 

g  =  0.01. 

while  the  initial  xj  and  X3  components  are  xj  =  X3  =  ^(1  -xf)/2.  The  increased 
accuracy  of  the  t  —*  00  phase  approximation  in  the  cartesian  solution  is  clearly  seen. 


5.1.2  Prolate  Gyrostats.  Observation  of  the  early  portion  of  spinup  in 
Figures  5.4  through  5.5  leads  us  to  believe  that  the  cylindrical  solution  is  better 
equipped  to  approximate  prolate  spinup  than  the  cartesian  solution.  Because  the 
cylindrical  solution  is  not  hampered  by  a  relatively  large  phase  error  in  the  early 
portion  of  the  spinup  maneuver,  we  expect  better  perform2uice  in  the  region  near 
the  instantaneous  separatrices. 

Figures  5.7  through  5.9  compare  the  numerical,  cylindrical  and  cartesian  solu¬ 
tions  for  spinup  of  a  prolate  gyrostat.  Here  we  see  the  large  influence  of  the  error 
made  by  the  xn  approximation  in  the  cartesian  solution.  In  essence,  the  cartesian 
solution  starts  out  with  the  wrong  information  and  the  solutions  keeps  getting  worse. 
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Figture  5.6.  Comparison  between  cartesian  and  cylindrical  equations  of  maximum 
error  detected  in  the  last  ten  percent  of  spin  up  as  a  function  of  initial 
conditions  for  an  oblate  gyrostat. 

The  cylindrical  solution,  however,  approximates  the  numerical  results  very  closely  in 
the  early  portion  of  the  spinup.  There  is  still  the  small  phase  error  that  occurs  at 
the  point  were  the  instantaneous  separatrix  crossing  really  does  occur.  Remember 
that  both  solutions  are  derived  from  linear  systems  that  approximate  the  nonlinear 
governing  equations.  The  total  absence  of  a  separatrix  in  either  approximate  solu¬ 
tion  is  the  cause  of  the  phase  error  that  develops  when  the  approximate  solutions 
attempt  to  model  the  governing  nonlinear  equations.  Figures  5.10  through  5.12  illus¬ 
trate  the  phase  error  in  the  cylindrical  solution  induced  by  inaccurate  modelling  of 
an  instantaneous  separatrix  crossing.  While  excellent  qualitative  results  are  evident 
in  Figures  5.7  through  5.9,  the  phase  error  at  the  instantaneous  separatrix  crossing 
propagates  the  error  in  the  solution  through  out  the  reminder  of  the  spinup. 

5.1.3  Transverse  Spinup.  Figures  5.13  through  5.15  compare  both  multiple 
scales  solutions  with  numerical  results  for  the  transverse  spinup  of  a  prolate  gyrostat. 
The  initial  conditions  are  slightly  offset  from  the  stable  flat  spin  equilibrium  point. 
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Figure  5.7. 


Figure  5.8. 
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Numerical  vs  cylindrical  and  cartesian  coordinate  approximate  solution 
of  the  xi  component  of  angular  momentum  of  a  prolate  gyrostat  with 
initial  conditions  Xi  =  0.4,  xj  =  0.64807,  X3  =  0.6^07  and  e  =  g  =  0.01 
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Numerical  vs  cylindrical  and  cartesian  coordinate  approximate  solution 
of  the  X3  component  of  angular  momentum  of  a  prolate  gyrostat  with 
initial  conditions  xi  =  0.4,  xj  =  0.64807,  X3  =  0.64807  and  e  =  g  =  0.01 
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Figure  5.9.  Numerical  vs  cylindrical  and  cartesian  coordinate  approximate  solution 
of  the  X3  component  of  angular  momentum  of  a  prolate  gyrostat  with 
initial  conditions  Xi  =  0.4,  xj  =  0.64807,  X3  =  0.64807  and  e  =  g  =  0.01 


Figure  5.10.  Percent  relative  error  of  the  Xi  component  of  angular  momentum  of 
a  prolate  gyrostat  using  cylindrical  coordinates  with  initial  conditions 
xi  =  0.7,  X2  =  0.4,  X3  =  0.59161  and  e  =  g  =  0.01 
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Figure  5.11. 


Figure  5.12. 
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Percent  relative  error  of  the  X}  component  of  angular  momentum  of 
a  prolate  gyrostat  using  cylindrical  coordinates  with  initial  conditions 
xi  =  0.7,  xj  =  0.4,  X3  =  0.59161  and  e  =s  g  =  0.01 
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Percent  relative  error  of  the  X3  component  of  angular  momentum  of 
a  prolate  gyrostat  using  cylindrical  coordinates  with  initial  conditions 
xi  =  0.7,  X2  =  0.4,  X3  =  0.59161  and  e  =  g  =  0.01 
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Figure  5.13.  Numerical  vs  cylindrical  and  cartesian  coordinate  approximate  solu¬ 
tion  of  the  xt  component  of  angular  momentum  of  a  prolate  gyrostat 
during  transverse  spinup. 

The  significant  increase  in  accuracy  during  the  early  portion  of  the  spinup  for  the 
cylindrical  solution  is  evident.  The  reason  for  the  improved  accuracy  is  the  same 
as  seen  during  prolate  spinup.  The  cylindrical  solution  is  not  hampered  by  a  larger 
error  early  in  the  spinup.  However,  the  essential  result  is  the  same  as  in  the  cartesian 
solution.  Extremely  small  torque  values  (which  implies  extremely  small  e  values)  are 
required  to  keep  the  trajectory  near  the  stable  flat  spin  equilibrium  point. 

5.8  Summary 

The  cylindrical  coordinate  multiple  scales  solution  developed  by  Lt  Col  Baker 
provides  excellent  results.  When  used  to  model  oblate  gyrostats,  the  approximate 
solution  is  much  more  accurate  in  the  early  portion  of  the  spinup  maneuver.  This  is 
due  to  the  fact  that  the  cylindrical  solution  is  not  dependent  on  the  xn  approxima¬ 
tion  that  the  cartesian  solution  has.  An  interesting  result,  however,  is  that  this  same 
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Figure  5.14. 


Figure  5.15. 
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Numerical  vs  cylindrical  and  cartesian  coordinate  approximate  solu¬ 
tion  of  the  xj  component  of  angular  momentum  of  a  prolate  gyrostat 
during  transverse  spinup. 
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Numerical  vs  cylindrical  and  cartesian  coordinate  approximate  solu¬ 
tion  of  the  X3  component  of  angular  momentum  of  a  prolate  gyrostat 
during  transverse  spinup. 
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Xii  approximation  provides  better  results  for  oblate  gyrostats  in  the  latter  portion 
of  the  spinup  maneuver. 


The  cylindrical  solution  does  have  superior  performance  when  modelling  pro¬ 
late  gyrostats.  While  the  solution  is  still  not  able  to  model  an  instantaneous  sep- 
aratrix  crossing  accurately,  it  does  provide  excellent  qualitative  information  for  the 
spinup  maneuver.  Further,  we  expect  that  spinup  which  begins  closer  to  the  prolate 
equilibrium  point  results  in  more  accurate  approximations. 

Finally,  both  the  cartesian  and  cylindrical  solutions  prove  to  be  inadequate 
when  modelling  transverse  spinup.  The  extremely  small  torque  values  required  ren¬ 
der  the  solution  impractical. 
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VI.  Conclusions  and  Recommendations 

6.1  Conclusions 

Two  approximate  solutions  for  the  spinup  of  a  near  axisymmetric  gyrostat  us¬ 
ing  the  method  of  multiple  scales  were  derived.  The  difference  in  transverse  inertias 
was  chosen  as  the  small  parameter.  This  selection  had  the  fortunate  property  of 
transforming  the  system  of  equations  that  described  the  spinup  from  a  nonlinear 
system  to  a  series  of  linear  systems.  The  first  solution,  derived  from  cartesian  co¬ 
ordinates,  incorporated  the  use  of  two  different  time  scales  to  capture  the  slowly 
varying  frequency  of  the  system.  To  derive  the  solutions,  an  approximation  had 
to  be  made  for  the  0{e)  term  of  the  xi  component  of  angular  momentum.  This 
approximation  induced  a  large  phase  error  in  the  transverse  components  during  the 
first  ten  percent  of  spinup  but  had  the  effect  of  reducing  phase  error  in  the  last  ten 
percent  of  spinup.  The  solution  was  compared  to  the  results  of  numerical  integration 
for  oblate  and  prolate  configurations.  Further,  the  solution  was  compared  for  the 
flat  spin  recovery  of  a  prolate  gyrostat.  Excellent  results  were  obtained  for  oblate 
configurations.  When  we  confined  our  error  analysis  to  the  last  ten  percent  of  the 
spinup  maneuver,  we  found  less  then  one  percent  error  in  all  three  components  of 
angular  momentum.  This  result  applied  to  gyrostats  with  a  one  percent  difference 
in  transverse  inertias  and  initial  conditions  in  the  northern  half  of  the  momentum 
sphere. 

The  results  for  prolate  configurations  were  hampered  by  the  extremely  non¬ 
linear  nature  of  prolate  spinup  trajectories.  Because  the  approximate  solution  was 
developed  from  a  series  of  linear  systems,  the  approximation  breaks  down  in  regions 
where  the  trajectory  is  rezJly  influenced  by  the  instantaneous  separatrices  associated 
with  prolate  spinup.  It  was  found,  however,  that  if  the  solution  of  the  numericad 
integration  at  a  point  in  time  after  the  last  bifurcation  on  the  momentum  sphere 
were  used,  reasonable  results  could  still  be  obtained  for  the  remainder  of  spinup. 
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Finally,  the  transverse  spinup  of  a  prolate  gyrostat  was  shown  to  be  completely  in¬ 
adequate.  The  trajectory  for  the  spinup  began  near  a  stable  flat  spin  center  and 
almost  immediately  crossed  2m  instantaneous  separatrix  due  to  an  excessively  large 
torque  value.  Reducing  the  torque  value  to  an  acceptably  small  value  would  require 
a  truly  axisymmetric  gyrostat. 

A  second  multiple  scales  solution  was  derived  by  Lt  Col  William  P.  Baker,  As¬ 
sociate  Professor  of  Mathematics,  United  States  Air  Force  Institute  of  Technology. 
Lt  Col  Baker  surmised  that  by  transforming  the  governing  equations  to  cylindri¬ 
cal  coordinates,  a  solution  could  be  derived  that  did  not  have  such  a  large  phase 
error  early  in  the  spinup  maneuver.  In  its  final  form,  the  cylindrical  solution  did 
not  have  an  explicit  slow  time  dependence,  but  did  have  the  slow  time  dependence 
show  up  implicitly  in  the  phase  equation.  A  comparison  between  the  cartesian  and 
cylindrical  solutions  for  oblate  gyrostats  produced  two  interesting  findings.  In  the 
early  portion  of  the  spinup,  the  cylindrical  solution  proved  far  superior  to  the  carte¬ 
sian  solution.  The  reason  was  that  the  cylindrical  solution  was  not  hampered  by 
any  approximations  in  its  derivation.  However,  in  the  final  portion  of  the  spinup, 
the  cartesian  solutions  phase  correction  at  t  =  oo  proved  more  accurate  than  the 
cylindrical  solution. 

The  cylindrical  solution  was  much  better  at  modelling  prolate  spinup  than 
the  cartesian  solution.  By  not  having  the  initial  phase  error  induced  by  the  xn 
approximation,  the  cylindrical  solution  was  better  at  modelling  the  trajectory  as  the 
instantaneous  separatrices  approached.  However,  the  solution  was  still  not  able  to 
accurately  model  the  instamtaneous  separatrix  crossing. 

6.2  Recommendations 

There  are  three  areas  in  which  further  study  is  warranted. 

1).  The  approximate  solutions  developed  in  this  thesis  were  only  derived  to 
approximately  0{e).  The  term  “approximately”  is  used  because  we  ignored  some  of 
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the  constants  of  integration  in  both  the  cartesian  2nid  cylindriceJ  solutions.  Addi- 
tionadly,  the  0{e)  expansion  limits  the  amount  of  torque  that  might  be  applied  to 
the  system.  If  more  accurate,  or  much  smaller  torque  values  are  desired,  the  sys¬ 
tem  of  equations  must  be  extended  to  0(e^).  This  will  require  either  an  analytic  or 
approximate  solution  to  the  following  integraJs. 


/  S.M 
jc.(c) 


COS  ada 
sin  ada 


2) .  For  prolate  configurations,  the  actual  trajectories  lie  in  regions  of  the 
momentum  sphere  influenced  by  instantaneous  separatrices.  If  the  initial  conditions 
start  fairly  close  to  the  prolate  equilibrium  point,  the  trajectory  encounters  smaller 
phase  error  at  the  separatrix  crossing.  A  matched  asymptotic  expamsion  might  be 
able  to  correct  the  approximate  solution  in  this  region. 

3) .  The  approximate  solution  derived  in  this  thesis  is  for  a  near  axisymmetric 
gyrostat.  Analysis  of  the  solutions  indicated  that  reatsonable  results  could  be  achieved 
for  a  gyrostat  with  a  two  percent  difference  in  transverse  inertias.  In  order  to  expand 
the  range  of  asymmetry  further,  am  alternate  technique  will  have  to  be  used  on  the 
nonlineau:  system  of  equations. 
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Appendix  A.  Straightforward  Expansion  Method  Solution  Terms 

The  following  terms  axe  to  be  used  in  Eqs.  3.32  and  3.33  for  the  approximate 
solution  using  the  straightforward  expansion  method.  Overbracketed  terms  are  sec¬ 
ular. 
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■*■  *3*/f^2aflc08^y^C08(/) 

~  *3*y|x5x3Cos(/)sm^y^ 
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Appendix  B.  Fresnel  Function  Integrals 


The  majority  of  equations  derived  in  this  thesis  were  obtuned  using  the  sym¬ 
bolic  manipulation  program,  Mathematica  2.0  for  SPARC  by  Wolfram  Research  Inc. 
(15).  This  program  was  able  to  solve  some  of  the  integrals  that  involved  the  Fresnel 
Sine  and  FVesnel  Cosine  functions.  However,  in  the  cylindrical  coordinate  derivation, 
the  program  was  not  able  to  solve  many  of  the  equations  and  so  they  were  solved 
by  Lt  Col  William  P.  Baker,  Associate  Professor  of  Mathematics,  United  States 
Air  Force  Institute  of  Technology.  These  derivations  are  included  here  in  full  as  a 
reference. 

We  begin  with  the  basic  definition  of  FVesnel  Sine  (S)  and  Fresnel  Cosine  (C) 
functions  (1:300) 


C{x)=  fo  coa 
S(x)=  /oSin(x4)df 


1)  The  first  integral  is  a  slightly  more  complex  cosine  function. 


g<T^  -I- 


2ka^d<r  =  j  cos  ^ 


da 


Using  a  simple  trigonometric  identity 


cos - 1-  sm 

9 


We  now  employ  a  change  of  variables.  Let 


V  = 

dv  = 


B-1 


The  restilting  integral  is 


which  has  the  solution 


C(i/)c08^ +S(i/)sin^ 


Using  a  similar  process  we  find  that 


it* 

S(j/)  cos - C{u)  sin  — 

9  9  J 


y/^iat+k) 


(B.l) 


(B.2) 


2)  The  second  integral  is  the  product  of  a  Fresnel  cosine  and  a  cosine  function. 

f  C^o(/(<^))  cos /(<y)d<r  (B.3) 

Jo 

As  in  Chapter  5,  we  define  the  following  terms 

Co(/)=  /  cos  fdt 
Jo 

•S’o(/)  =  /  sin/d< 

Jo 

where 

f  —  gt^  +  ikt  +  $0  (®*^) 

Note: 


=  sin/(t) 
=  cos  fit) 


B-2 


To  solve  the  integral  we  employ  the  integration  by  parts  technique,  let 


«  =  cum) 

du  =  cos  f{t)dt 


L 


udu  =  —u^ 
0  2 


Using  a  similar  process  we  find 


/„ 


‘  = i  [c.  (/(())r 
0  2 


1 


SUm))s\nf{<T)da  =  -[SUm)Y 

0  £ 


3)  The  third  integral  is  a  complicated  Fresnel  cosine  function. 


jf  +  2ka^  da 


Using  the  definition  of  a  Fresnel  cosine  we  have 


-a:  cos  +  2ka^dads 


Changing  the  order  of  integration  results  in 


=  X/;  cos  +  2ko^  dsda 
=  j  {t  — a)  cos  +  2ka^da 


Now  we  add  and  subtract  the  same  term 


(B.5) 


(B.6) 


=  t  y  cos  ^ga^  +  2ka^da  ~  j  +  2fc)  cos  (jga^  +  2ka^dc 


+  —  y  cos  (^ga^  +  2ka^da 


B-3 


The  resulting  solution  is 


=  (gt^  +  2kij  —  i  sin  (jgt^  +  2kij  +  -Co  (gt^  +  2kij  (B-7) 


Using  a  similar  process  we  find 


J  So  {^(T^  +  2k<r^  da  =  tSo  (^gt^  +  2kt^  +  ^  [cos  +  21:f)  —  l|  +  -  So  (gt^  +  2kt^ 

(B.8) 

4)  The  fourth  integral  involves  the  product  of  the  function  and  the  variable  of  inte¬ 
gration. 

j  a  cos  (f{a))da  =  j  a  cos  (ga^  -f  2ka  -f-  Bo^da  (B.9) 

Agsdn,  we  add  and  subtract  the  same  term 

1 


=  —  y  {2ga  +  2A:)  cos  ^ga^  -I-  2ka  -I-  Oojda 
—  —  J  cos  -h  2ka  -f  Bo)da 


The  solution  to  the  integral  is 


=  i[»m/(i)-sm<lo]-|a  (/(!)) 


Using  a  similar  process  we  find 


f  as\n{f{a))da  =  [cos /(f)  -  cos^o]  -  -S<,  (f(t)) 
Jo  £g  g 


5)  The  fifth  integrsJ  is  the  cosine  function  squared. 


/t  1  r* 

/  cos^  f  da  =  -  /  [1  -f-  cos  2/]  da 
Jo  2  Jo 


=  I +  5^.(2/) 


(B.IO) 


(B.ll) 


(8.12) 


B-4 


6)  The  sixth  integral  is  the  product  of  the  variable  of  integration,  a  Fresnel  sine 
function  and  a  cosine  fimction. 


l\{2g<T  2k)  S,{f)  cos  f]  da 

Jo 

We  employ  the  integration  by  parts  technique.  Let 

ti  =  SM) 

dv  =  {2ga  +  21:)  da 

This  results  in 

=  So{f)  sin  /lo  -  /  sin^  fda 
Jo 

Using  a  trigonometric  relation,  we  have  the  solution 

=  5,(/)8m/-i|  +  ic.(2/) 


B-5 


(B.13) 
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